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1.0 SUMMARY 


This document describes the analysis and use of program L221 (TEV156), a digital 
computer program to calculate steady state solutions for linear second-order differential 
equations caused by sinusoidal forcing functions. The program calculates generalized 
coordinate frequency responses, load frequency responses, and load power spectral 
density (PSD) parameters using random harmonic analysis techniques. The particular 
field of application of the program is the analysis of airplane response and loads caused 
by continuous air turbulence, both vertical and lateral gusts. 

The equations of motion and load equations of the vehicle must be formulated outside 
L221 (TEV156) and read as program input data from either cards or magnetic files. The 
frequency response functions for the generalized coordinates and loads (shears, 
moments, accelerations, etc.) are solved at a large number of discrete frequencies. An 
empirical representation of the atmospheric turbulence PSD is utilized together with 
the load frequency response functions to obtain the load PSD pa'rameters A and No (root 
mean square (RMS) load/root mean square gust velocity and the number of zero 
crossings with positive slope per unit distance, respectively). 

A number of options are available within the program to produce flexibility as to the 
type of solution and results desired. Some of these features are: 

• Equations of motion and load equation coefficient matrices may be either frequency 
dependent (nonconstant coefficients) or frequency independent (constant 
coefficients) 

• User selection of gust spectra description 

• User-specified scale of turbulence 

• Multiple forcing function spectra description 

• Static-elastic solution for quasi-steady path 

• Deletion of user-specified degrees of freedom 

• User-specified structural damping factor for each degree of freedom 

• User-specified starting point for RMS integration 


• User-specified Wagner and Kiissner indicial lift growth function coefficients 

• Gradual gust penetration 

• Load correlation 

• User input or internal generation of the array of solution frequencies 

The following results can be saved on magnetic tape for automatic plotting with 
user-supplied interface and plotting programs: 

• Frequency response functions for the generalized coordinates* 

• Frequency response functions for the loads* 

• Load power spectra functions* 

• Real and imaginary components of the generalized coordinate frequency responses 

• Real and imaginary components of the load frequency responses 

• RMS loads, A* 

*These items can also be automatically plotted using the COMpSO plotter. 

Limitations imposed on the equations are: 

• All roots must have positive damping. 

Note: A solution will occur even though the damping is negative; however, it will 
be in error as no steady state solution exists. 

• The matrix which is the sum of the zero-, first-, and second-order coefficients must 
be nonsingular at all solution frequencies. 
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2.0 INTRODUCTION 


The computer program L22KTEV156) may be used as either a stand-alone program or 
as a module of a program system called DYLOFLEX which was developed for NASA 
under contract NASI- 13918 (ref. 1). Because of the DYLOFLEX contract requirements 
developed in reference 2, a program was needed to calculate power spectral density 
(PSD) gust load parameters for equations of motion and load equations which have 
frequency-dependent coefficients. An existing program^ that calculated PSD gust load 
parameters for equations with constant coefficients was modified according to the 
DYLOFLEX specifications,^ to optionally use nonconstant coefficient matrices for the 
equations of motion and load equations. 


^Clemmons, R. E.: A Power Spectral Digital Computer Program to Determine Dynamic Loads due 
to Random Gusts-PSDSYS (TEV 1 56)-U8ers Guide, BCS-G0235, June 1973. 

2 

Clemmons, R. E.: Programming Specifications for Modules of the Dynamic Loads System to 
Interface Vlith FLEXSTAB. NASA contract NASl-13918, BCS-G0701, September 1975. 
(Internal Documents) 
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3.0 SYMBOLS AND ABBREVIATIONS 


The following list contains items that appear in this document except section 6.3 (card 
input). 


Engineering 

notation Definition 


{A} 

[A].[B] 

[AIC] 

ai.bl,ci,di 

a2.b2>C2>d2 

c 

C 

C 

{C2}.[C2] 

{C2}.[C2] 

{Cal.LCs] 

{C3}.[C3] 

{CS0}.{CSl},{CS2} 

{Fa} 

{f^} or . I 

{FREQM} 

f(s) 


Gust response factor = o-f/o-g, 
response/unit velocity. 

Real and complex coefficients of the 
AIC matrix. 

Aerodynamic influence coefficient 
matrix. 

Wagner indicial lift growth 
coefficients. 

Kiissner indicial lift growth 
coefficients. 

Reference chord length, units of length 

Scalar constant on the equations of 
motion forcing function when using 
gradual penetration. 

Scalar constant on the load equation 
forcing function when using gradual 
penetration. 

Equations of motion forcing function 
coefficient matrix. 

Load equations forcing function 
coefficient matrix. 

Equations of motion convoluted forcing 
function coefficient matrix. 

Load equations convoluted forcing 
function coefficient matrix. 

Special feedback coefficients. 
Aerodynamic force. 

Streamwise distance from the point 
first encountering gust to the points 
encountering the gust later, units of 
length 

Array of frequencies at which the 
frequency-dependent input matrices 
are defined. 

Laplace transforni of the Kiissner 
indical lift growth function. 
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Fydamp 

Special feedback transfer function. 

G 

Feedback gain coefficient of the 
feedback transfer function. 

{ksdI 

Structural damping coefficient matrix. 

g(s) 

Laplace transform of the Wagner 
indicial lift growth function. 

i 

Vrr 

L 

Scale of turbulence, units of length. 

[M] 

Lumped inertial data. 

[M,].[M2].[M3] 

Equations of motion coefficient 
matrices for generalized coordinate 
displacement, rate, and acceleration, 
respectively. 

[M4].[M5].[M6], 

Equations of motion coefficient 
matrices for convoluted generalized 
coordinate displacement, rate, and 
acceleration, respectively. 

[Mi].[M2].[Si3] 

Load equations coefficient matrices for 
generalized coordinate displacement, 
rate, and acceleration, respectively. 

[M4],[M5],[M6] 

Load equations coefficient matrices for 
convoluted generalized coordinate 
displacement, rate, and acceleration, 
respectively. 

{No} 

Number of zero crossings with positive 
slope per distance. 

{P} 

Coefficients of the special feedback 
transfer function. P(1),P(2),P(3) and 
P(4),P(5),P(6),P(7) are the second-order 
functional coefficients for the 
numerator and denominator, 
respectively. 

q 

Dynamic pressure, fore/unit area 

{Q} or {q) 

Generalized coordinates of the 


equations. 

QMAGl 

Magnitude of the load solution at each 
solution frequency. 

QMAGr 

Magnitude of the generalized 
coordinates at each solution frequency. 

{R} 

Coefficients of the special feedback 
transfer function. R(1),R(2),R(3) and 

\ 

R(4),R(5),R(6),R(7) are the zero-order 
functional coefficients for the 

1 

numerator and denominator, 
respectively. 
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[Si].[S2],[S3] 

[S4],[S5].[S6] 

[Si],[S2],[S3] 


[S4].[S5].[S6], 

[SPEC] 

{SUMC} 

{SUMC} 

[SUMM] 

[SUMM] 

{T} 


{u} 


Vt 

{W} 

{Wgi} 

T(ift) 

Z 


Laplace operator = io). 

Equations of motion special feedback 
coefficient matrices for generalized 
coordinate displacement, rate, and 
acceleration, respectively. 

Equations of motion special feedback 
coefficient matrices for convoluted 
generalized coordinate displacement, 
rate, and acceleration, respectively. 

Load equations special feedback 
coefficients matrices for generalized 
coordinate displacement, rate, and 
acceleration, respectively. 

Load equations special feedback 
coefficient matrices for convoluted 
generalized coordinate displacement, 
rate, and acceleration, respectively. 

Load output spectrum at each solution 
frequency. 

Equations of motion forcing function 
summation matrix. 

Load solutions at each solution 
frequency. 

Generalized coordinate summation 
matrix. 

Load equations summation matrix for 
generalized coordinates. 

Coefficients of the special 

feedback transfer function. T(l) 

and T(2),T(3),T(4),T(5),T(6),T(7) 

are first-order functional coefficients 

for the numerator and denominator, 

respectively. 

Coefficients of the special feedback 
transfer function. U(1),U(2),U(3), and 
U(4),U(5),U(6),U(7) are first-order 
function coefficients for the numerator 
and denominator, respectively. 

True airspeed velocity, length/sec. 
Normalwash. 

Normalwash caused by the gust from 
point i. 

Complex frequency response function. 
Scale factor on forcing function. 

Wagner indicial lift growth 
coefficients. 





♦(t) 

<Mt) 

iftg 


Kiissner indicial lift growth 
coefficients. 

Gust angle, radians. 

Wagner function. 

Kiissner function. 

Mode deflection 
Mode slope. 

Gust mode shape. 

Root mean square (RMS) of response. 
Root mean square of gust velocity with 
zero mean. 


<b(n) 

n 

tu 

£ 

[^1 

m 

i 

®*^llm 


K 


V 


u 

[] 

n 



♦ 

A 


Gust spectrum. 

Spatial frequency = (o/V, radian/unil len^h 
Frequency, rad/s. 

Laplace transform. 

Equations of motion gradual 
penetration forcing function matrix. 

Load equations gradual penetration 
forcing function matrix. 

Structural damping factor. 

Rudder angle. 

Rudder limit for special feedback 
transfer function. 

Constant for the exponential input 
spectrum. 

Exponential power constant for the 
exponential input spectrum. 

Row matrix. 

Rectangular matrix. 

Column matrix. 

Diagonal matrix. 

Transpose of the matrix. 

Indicates convolution. 

Real part of complex number. 



4.0 ENGINEERING AND MATHEMATICAL DESCRIPTION 


This program is capable of calculating, steady state solutions for any linear second-order 
differential equations caused by sinusoidal forcing functions. Statistical characteristics 
of loads for the forcing function spectral shape may also be calculated using random 
harmonic analysis techniques. 

The following engineering description is for the specific field of application: the analysis 
of airplane response and loads due to continuous random atmospheric turbulence. 

4.1 EQUATIONS OF MOTION 


In general, the equations of motion consist essentially of three parts: 

1. Structural (generalized stiffness, inertia, and damping) 

2. Aerodynamic (generalized response and excitation air forces) 

3. Stability augmentation system (SAS) (equations relating system transfer functions 
and sensor signals) 

The equations of motion can be represented in the form: 

(Structural + SAS) + (Response Aerodynamics) = (Gust Excitation Aerodynamics) 


([Mi]{q} + [M 2 ]{q} + [Majfq}) + ([M4]{q}*«I> + [M5]{q}*4>) 

= -Z ({C2}ag + {Cajag*!/;) 

where: 

M’s, C’s = Appropriate matrix coefficients 

q’s = Generalized coordinates including SAS degrees of freedom 

«g = Gust angle 

4> = Wagner function (equal one for no lift growth) 

ifr = Kiissner function (equal one for no lift growth) 

* Indicates convolution 

The composition of [Mj], [M 2], and [M3] can be visually partitioned into the 
structural and SAS submatrices 

Mi I 

^structi 



where: 


^^Istruct^ 

^^2struct^ 

^^^Struct^ 


= Generalized structural stiffness 
= Generalized structural damping 
== Generalized structural inertia 


The generalized structural inertia structural 

stiffness “ ["^] [^]'*'[^] [^] usually formed in vibration programs 

which require only geometric, structural stiffness, and mass data. 


Generally, structural damping is assumed to be proportional to displacement but in 
phase with velocity. It can be written in the present notation as: 

»[BSD][Migjpm.j]{q} (3) 

where ggo = 2 (structural damping factor, on each degree of freedom). 

This representation is generally associated with harmonic motion and is built into this 
program as an option. However, if an equivalent viscous damping representation is 
desired, the structural damping can be transformed into an equivalent viscous damping 
representation given as, 


[kd][m,;^^,] "*{4} (4) 

This representation is sometimes used with harmonic motion. However, it should be 
noted that when this method is used, the structural damping factor will vary with 
frequency and have a value of gsD^^ natural frequency of each mode. 


The coefHcients 

not a function of the aerodynamic representation. Therefore, they are placed in the Mi 
through M3 matrices. 


ISAS. 


'Ji=l 3 the deflnition of the SAS equations and are 


Because this is a general-purpose program, the M4 and M5 matrices usually consist of 
the generalized response aerodynamic forces, with that portion of M4 and M5 
corresponding to the SAS submatrices in Mi, M2, and M3 equal to zero. The generalized 
response aerodynamic matrices are formed by the expression; 


where: 

{Fa} = Aerodynamic force 
[<^] = Mode deflection 
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In general, the aerodynamic force is given as 

{Fa} = q[AIC]{W} (5) 

where: 

{W} = Normalwash vector i 

- ^ (*){q} 

q = Dynamic pressure 

[<^'] = Mode slope 

[AIC] = [a] + i[B]; [B] = 0 if (reduced frequency) k = 0 ^ 

Note: M4 and M5 are only valid for the specific w corresponding to the reduced 

frequency, k, that was used to calculate the imaginary AIC matrix [B]. 

The generalized excitation matrix {C3}, which is generally caused by gust, may remain 
complex and similarly may be written: 

{Cgjag = q[^]'^[AIC]{Wg} | (6) 


where (Wg) = {<^g}e*^* for harmonic motion, 
also 


then 


so define 


w« = e“ (‘- V ) 

iwt V 
= e e * 


■(*•1 


Note: If Jtj = 0 (no gust penetration), then {<^g} = {l}. 


(7) 


4.2 LOAD EQUATIONS 

The load equations follow the same format as the equations of motion. 
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{Load} = [Mi]{q| + [M2](q} + [MsJjq} 

+ [^]{q}** + + |z|({C2}«g + {CslV'*') 

where: 

[Mi],[M 2],[513] = Load matrix coefficients of the generalized coordinate displacement, 
rate, and acceleration, respectively 

= Load matrix coefficients of the generalized coordinate rate and 
acceleration convoluted with the Wagner function 

= Load matrix coefficient of the excitation function (usually gust 
angle ag) 

= Load matrix coefficient of the excitation function convoluted with 
the Kiissner function 

= Wagner function 

= Kiissner function 

- Generalized coordinate displacement, velocity, and acceleration 
responses, respectively 

Indicates convolution 

Relating these matrices in a physical sense, [M|], [MoL and [M3] are usually 
associated with the load resulting from structural response; [M4] and [^5] are usually 
associated with the load resulting from aerodynamic response; {^3} is usually 
associated with the load resulting from the gust excitation force. 

4.3 SOLUTION TECHNIQUE 

The steady state solution for the equations of motion and load equations caused by 
sinusoidal forcing functions is obtained by using the Laplace transform technique. 
Taking the Laplace transform of the equations of motion and load equations results in 
the form: 

Equations of Motion 

[Ml] + + s[M2] + b 2[M3] + g(s)[M4] + 8g(8)[M5] [ir(q)} 


[M4],[^5] 

{^2} 

{C3} 

<I> 

q.q.q. 


= *Z [{C 2 } + f(8){C3}]jC(ag) 


( 10 ) 



where: 


g(s) 

f(s) 


= sjC |$] = - 

= s£ [<//] =32- 


bjS 

CjS 

d]S 

S + Q!j 

s+ 

s + 7i 

b2S 

C2S 

d2S 

8 + ^2 

s + ^2 

's + 72 


KSD == Structural damping factor 

Z = Scaling factor 

s = Laplace operator = io> 

and using the relationship that the Laplace transform of the convolution theorem gives 
Duhamel’s formula: 




ft 

=fo' 


=X /f(t - x)q(x)dx 




J 


( 11 ) 


Load Equations 


Load 


= {[Mi] + s[M 2 ] + s^[M 3 ] + g(s)[M 4 ] + sg( 8 )[M 5 ]|{£(q)} 


+ |Z|{{C2} + f(s){C3ll£(ag) 


( 12 ) 


Thus, generalized coordinated and load frequency response functions are easily obtained 
by solving equations (10) and (12). 

The matrix coefficients for the equations of motion and load equations may be either 
independent of frequency or frequency dependent. However, for the frequency-dependent 
matrices, a linear interpolation of the matrix coefficients is used to form additional 
matrices at frequencies intermediate to those formed external to this program. This 
interpolation process is necessary so that the equations of motion and load equations 
may be evaluated at a sufficient number of frequencies to adequately define the 
frequency response function, while generating as few input matrices as necessary. The 
selection of frequencies is dependent upon the damping in the system, the frequency 
range over which the contribution of the responses significantly affect the final RMS 
loads, the number of modes used, and the frequency separation between modes. For 
typically low-damped dynamic models, this requires the solution of the equations of 
motion at a large number of frequencies (on the order of 200 to 250). The selection of 
frequencies is delegated to the engineering user, and the number must be kept to a 
minimum to reduce costs, yet must be sufficient to adequately define the frequency 
response functions. This interpolation is required for only the response and gust 
excitation aerodynamic matrices since the structural matrices and the SAS equations 
remain constant with frequency. 
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Power spectral density (PSD) load parameters, root mean square (RMS) load/root mean 
square gust velocity 0^), and number of zero crossings (Nq) can be calculated using the 
load frequency response functions. The equations for X and Nq are 


and 



|T(ift)|^ d)(n)da 


IP 


1l/2 


(13) 


No 


1 

2irA 


00 

/ ft*|T(ift)|^«D(n)dft 
0 


1l/2 


j 


where: 

T(ift) = Load frequency response 
<l>(fl) = Gust spectrum 


(14) 


In practice, the upper limit (infinity) of the integral is replaced by the finite value of 
the last discrete frequency used in the solution of the transfer function. 
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5.0 PROGRAM STRUCTURE AND DESCRIPTION 


Program L221 (TEV156) has been constructed as an overlay system consisting of a main 
overlay and three primary overlays: 

Main overlay (L22 1,0,0) L221 

Primary overlay (L22 1,1,0) FINDRMS 

Primary overlay (L221,2,0) ^SORTQLS 

Primary overlay (L22 1,3,0) PLOTQLS 

The main overlay L221 simply reads cards to determine which primary overlays are to 
be called, calls the proper primary overlays into execution, and aids in communication 
between overlays through labeled common blocks. 

FINDRMS performs the actual calculations of L221 (TEV156). 

SORTQLS sorts the results of FINDRMS and prepares a magnetic file of data vectors 
suitable for plotting. 

PLOTQLS sorts the results of FINDRMS and prepares a magnetic file of plot 
instructions for the COMpSO plotter. 

A schematic of the L221 (TEV156) overlay structure is presented in figure 1. A 
description of each overlay follows. 

5.1 PROGRAM L221 

L221 is the program name of the (L221,0,0) overlay. This overlay repeatedly reads cards 
to determine the type of execution to be performed and calls primary overlays to 
perform the necessary operations. 

The following is a list of keywords used to control the execution of the system’s primary 
overlays. 

$FREQuency Indicates that the input data following this card is for use by 

L221 (TEV156). 

$TITLe Is printed and ignored. 

$FINDrms Causes the execution of program FINDRMS, (L221,l,0), which 

performs all of the computations of L221 (TEV156). 

$SORTqls Causes the execution of program SORTQLS, (L221,2,0), which 

sorts the results of FINDRMS and prepares a magnetic file of 
data vectors suitable for plotting. 
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A's and Zero Crossings 


All magnetic files except "TAPE99" are written in READTP/WRTETP formats. "TAPE99" 
is a file written by PLOTQLS which is suitable for plotting on the COMpSO. 


Figure 1.—L221 (TEV156) Overlay Structure and Input/Output Fites 
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$PLOTqls 


$CORE 

$CHECk 


$QUIT 


Causes the execution of program PLOTQI^, (L 221 , 3 , 0 ), which 
sorts the results of FINDRMS and prepares a magnetic file of 
plot instructions for the COMpSO plotter. 

Causes the maximum field length limit to be changed. The field 
length limit default value is 200000 octal. 

Triggers a special option to aid a programmer testing fatal 
error logic paths within L221 (TEV 156 ). The option causes the 
clearing of the fatal error flag and allows execution to continue. 
This option is not intended for use during normal program 
executions. 

Causes L221 (TEV 156 ) to terminate execution. 


5.2 PROGRAM FINDRMS 

Program FINDRMS, overlay (L221,l,0), performs all of the computations of L221 
(TEV 156 ). Given an input power spectrum 4 >(ft) over a range of frequencies 
(ftl, 02, . . OnfREQ)> matrices of equations of motion, load equations, and excitation 
forces, FINDRMS will find: 

{Q} Response function of the generalized coordinates at each solution frequency 

{SUMC} Response function of the load at each solution frequency 

and integrate over the range of frequencies to find: 

{a} RMS load values 

{Nq} Number of zero crossings per unit length 

The FINDRMS operations are divided into the following 10 steps to obtain the solutions 
of the generalized coordinates. 

1 . Read card input of options and constants. 

2 . Read card or magnetic file input of equations of motion [Mi], [M2], [Mg] and 

[Si], JS2], [S0] for feedback plus the excitation forces {C2} and {C3} or {f^} 

and [<^] for gust penetration. This includes {FREQM} (frequencies associated with 
k- value aerodynamics) if NKVAL > 0 . 

Note: [Mq], [Si], and {C2} are not used when using the DYLOFLEX system. 

3 . For each of the NFREQ frequencies, solve for {Q} the generalized coordinates, save 
{Q} on the scratch file, and optionally print {Q}. The generalized coordinates are 
found by solving the following equation. 

[SUMM] {Q} = {SUMC} ( 15 ) 
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where: 


[SUMM] = [M,] + i[gSD][Mi] + s[M2] + 8^[M3] 

+ g(s)[M4] + sg(s)[M5] + 8^g(s)[M6] ( 16 ) 

where; 

s = ico, Laplace transform variable 


g(s) = Laplace transform of the Wagner function 
b j s c I s d j s 
s + Oij s + |3j s + 7| 


For a solution with feedback, the following terms are also added to [SUMM]. 

Fydamp |[Si] + s[S2] + s^[S3] + g(s)[S4] + sg(s)[S5] + s^g(s)[S6| ( 17 ) 

where: 


Fydamp = 


G(Pxs^ + U^s + R|)(p2S^ + U2S R2)(p38^ + U3S + R3)(T]S + 1) 
(P4S^ + U4S + R4),..(P78^ + U78 + R7)(T28 + 1)...(T78 + 1 ) 


{SUMC} =-zto{C3} + {C2}} 


( 18 ) 


where: 


f(s) = Laplace transform of the Kiissner function 

^2^ 


^2^ 


^2^ 


^2 s + 0f2 ^ ^ "^2 


When gradual penetration is requested, {C3} is frequency dependent and will 
be calculated as 


{C3} = C <[$]cos{ft{fjt}} - i[$]sin{n{f^}}l 


( 19 ) 


Note; 


Optionally, the input matrices [M4], [M5], [Mg], and {C3} or [?i are 
frequency dependent. Interpolation is required to define those matrices at 
each frequency (NFREQ), before calculating [SUMM], {SUMC}, and {Q}. 

Optionally, the user may specify multiple forcing functions causing the 
matrices {C3}, {SUMC}, and {Q| to be rectangular. (See appendix C for a 
discussion of the multiple forcing function option.) 

Optionally, the user may specify certain degrees of freedom to be deleted. 
This is accomplished by zeroing the appropriate degrees of freedom in the 
[SUMM] and {SUMC} matrices. 
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4 . Modify the equations of motion for a static-elastic solution and repeat step 3 . 
Chosen columns of the matrices [M2], [M3], [M5], and [M0] are set equal to zero 
for the static-elastic solution; optional. 

The loads are solved by the following steps. 

5 . Read card_or magnetic file input of load equations [Mi 2 i [M2L ...,[M0] and [Sj], 
[S2], [Sg] for feedback plus the matrices {C2} and {C3} or [<^»]. 

Note: [Mg], [^j], and {C2} are not used when using the DYLOFLEX system. 

6. For each frequency read {Q} (fro m step 3 ); calculate {SUMC}, the input spectrum 
0(n), and the output spectrum {SPEC}; optiona lly wri te data on magnetic files 
"'IPLTPE” and ”IRTAPE”; optionally print {SUMC}, and keep a running 
integration over the frequencies of {A} and {Nq}. 


where: 


{SUMC} = [SUMM]{Q} +|Z|{f(s){C3} + {C2} 


[SUMM] = [Ml] + s[M2] + s^[M3] + g(s)[M4] + sg(s)[Ms] + sVs)[Mg] 

For a solution with feedback, the following terms are also added to [SUMM]. 

Fydamp [§i] + s[S2] + s^[S3] + g(s)[S4] + sg(s)[S5] + sVs)[S6]] ( 21 ) 

When gradual penetration is requested, {C3} is frequency dependent and will be 
calculated as: 


{C3} = c|[$]cos|n{fji}|- i[ 5 ]sin|n{f£^^ 
The integration for {A} and {Nq} at one load, ^ , is shown as: 


' NFREQ / SPECgi + SPECgj.i 




/ , \ NFREQ / 

(^) V ifi ( 


NFREQ / sPECa £2;^ + SPECji., nj.,; 


) (S2i-ni.i; 


where SPECj^l is the output spectrum calculated for load i at frequency i and is 
found by: 

SPECjii = 4 >(Cli)* 1 SUMCju I ^ (25) 

where |sUMCjii| is the absolute value of the complex number SUMC 
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The input spectrum value <^(n) is found by the user’s choice of one of the three 
following equations: 


IT 


Dryden 


1 +3(Ln)^ 

[1 +(LS2)2]2 
1 + 8/3(1. 339LQ)2 
[1 +(1.339LSl)2]-y- 


Von Karman 


0(S2) = K (Siy^ Exponential 


(26) 

(27) 

(28) 


or by interpolation over a user-supplied tabular spectrum. 

7. Print {A} and {Nq} and optionally punch them on cards. 

8. Check the correlation between the different load responses requested by card input 
data; optional. 

Note: Optionally, the input matrices [M4], [M5], [M0], {C3} or [0] are 
frequency dependent. Interpolatio n is req uir ed to d efine these matrices at 
each frequency before calculating [SUMM], {SUMC}, etc., in step 6. 

Optionally ,_the user may s pecify multiple forcing functions causing the 
matrices {C3} and {SUMC} to be rectangular. (See appendix C for a 
discussion of the multiple forcing function option.) 

Optionally, the user may specify certain degrees of freedom to be deleted. 
This is accomplished by zeroing the appropriate degrees of freedom in the 
[SUMM] matrix. 

9. Modify the load equations for a static-elastic solution; set chosen columns of [M2], 
[M3], [M5] and [Mg] to zero; and repeat steps 6, 7, and 8; optional. 

10. For additional sets of load equations, repeat steps 5 through 9; optional. 

5.3 PROGRAM SORTQLS 

Program SORTQLS, overlay (L221,2,0), which sorts generalized coordinate and load 
responses, is called to process data written on the magnetic file *TPLTPE” by program 
FINDRMS and to write an output magnetic file "NEWTPE” which contains pairs of 
matrices, independent and dependent variables, to be plotted by a subsequent program. 

program. 

"IPLTPE” contains one file of data (see fig. 11) for each FINDRMS solution. Each file 
begins with the frequency array {FREQ}, a 1 x NFREQ matrix, and is followed by 
NFREQ groups of the three or four arrays below. (A) will be present only if 
IPLRMS = 1 on FINDRMS card 4.3. 
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6.0 COMPUTER PROGRAM USAGE 


The program was designed for use on the CDC 6600. The machine requirements to 
execute L221 (TEV156) are: 

Read control cards and card input data. 

Print standard output information, optional intermediate 
calculations, and diagnostic messages. 

All magnetic files not specifically defined as magnetic tapes are 
assumed to be disk files. 

Tape drive For permanent storage of data. Magnetic files are copied to and from 

magnetic tapes with control cards before and after program 
execution. 

Card punch Punch the final X values and number of zero crossings; optional. 

Program L221 (TEV156) is written in FORTRAN and may be compiled with either the 
RUN or FTN compiler. L221 (TEV156) may be executed oh either the KRONOS 2.1 or 
NOS operating system. 

6.1 CONTROL CARDS 

The following list is a typical set of control cards used to execute L221 (TEV156) using 
the absolute binaries from the program's master tape. 

Job Card 
Account Card 


REQUEST (MASTER, F=I,LB=KL,VSN=66XXXX) 
REWIND (MASTER) 

SKIPF (MASTER) 

COPYBF (MASTER, L221) 

RETURN (MASTER) 

• 

L221. 

EXIT. 



Card reader 
Printer 

Disk storage 
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DMP (0, field length) 
— End-of-record 


' Card Input 
Data 


— End-of-file 


The following list is a typical set of control cards used to execute L221 (TEV156) using 
the relocatable binaries from the program’s master tape. 

Job Card 
Account Card 


• 

REQUEST (MASTER.F=I,LB=KL,VSN=66XXXX) 
REWIND (MASTER) 

SKIPF (MASTER.2) 

COPYBF (MASTER, REL221) 

RETURN (MASTER) 

• 


LDSET (LIB=DYLIB,PRESET=INDEF) 
LOAD (REL221) 

NOGO. 

RETURN (REL221.DYLIB) 

L221. 


EXIT. 

DMP (0,field length) 
— End-of-record 




Retrieve the program 
from its master tape 


Prepare optional 
input files and 
i retrieve DYLIB, the 
DYLOFLEX alternate 
^ subroutine library 

j Load and execute 
1 L221 (TEV156) 


{ Save optional output 
data files 


^Card Input 
Data 


— End-of-file 


6.2 RESOURCE ESTIMATES 

The computer resources used (core requirements, tapes, printed output, time, etc.) are 
functions of the problem size and program options used. Table 1 contains examples of 
the resources used for a number of sample problems. 
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Table 1.— Resources Used by Sample Problems 



NDOF 

NFREQ 

NFORC 

NKVAL 

NLD 

CPU Seconds 
Disk Requests 
Disk Sectors 


Required Core; 
(Octal ) 

FINDRMS 

SORTQLS 

PLOTQLS 


SAMPLE PROBLEMS 

1 

2 

3 

4 

5 

7 

7 

7 

8 

7 

10 

10 

10 

210 

72 

0 

0 

0 

0 

6 

0 

0 

10 

20 

0 

6 

6,9* 

6 

8 

7 

2 

25 

2 

15 

10 

100 

450 

80 

310 

350 

830 

3500 

350 

1300 

1800 

800 

2000 

200 

2400 

3300 

122000 

122000 

122000 

123000 

125000 





124000 


43000 





66000 





The problem has two load sets 








F&eld Length 


The field length required by L221 (TEV156) is dependent on the problem size and the 
program module(s) used. Core must be requested for the largest module to be run; either 
FINDRMS, SORTQLS, or PLOTQLS. The default majnmum core size established in 
L221 (TEV156) is 200000 octal. 

Field length = FLp + FLm 

where FLp is the program length, and FLm is the matrix storage required. 

The module(s) have the following octal values of FLp: 

FINDRMS 126000 

PLOTQLS 66000 

SORTQLS 60000 

For FINDRMS, the computed value of FLm is the larger of FLq or FL £. 

FLq = Core required for equations of motion 
= NFREQ 

+ NDQF(8 (NFORC) + 3) 

+ NDOF^(NMi + 2) 

+ NDOF; if static elastic 
+ NDOFD; if deletion of freedoms 
-f- NDOF(2 (NPAN» + NPAN; if gradual penetration 
+ NDOF^(NSi) + 3(NDOF); if feedback 
+ NDOF^(4 (NFORC) + 6) + NKVAL; if NKVAL > I 
+ NDOF(4 (NPAN)); if NKVAL > 1 and gradual penetration 
FL I = Core required for load equations 
= NFREQ + 3(NDOF) 

+ 3(NLD) 

+ NLD(NDOF)(NMi + 2) 
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+ NLD(6 (NFORO) 


+ NFORC(NFORC + I) 

+ NDOF; if static elastic 

+ NDOFD; if deletion of freedoms 

+ 32(ICORR); if correlation 

+ NPAN(2 (NLD) + 1); if gradual penetration 

+ NLD(NDOF)lNSi) + 3(NDOF); if feedback 

+ NLDC6 (NDOF) + 2(NFORC) + NKVAL); if NKVAL > 1 

+ NLD(4 (NPAN)); if NKVAL > 1 and gradual penetration 

For SORTQLS, the computed value of FLm is: 

FLm = NFREQ*(NR + 1) + NFREQ*NR; if pairing is to take place 

where NR is the larger of NDOF or NLD. 

For PLOTQLS, the computed value of FLm is the same as for SORTQLS. 

Note: FLm for FINDRMS, SORTQLS, and PLOTQLS must be converted to an octal 
base number before being added to FLp to find the required field length. 

Time Estimate 

The following factors will be used to calculate the CPU time in seconds for L221 

(TEV156), FINDRMS. SORTQLS. and PLOTQLS. 

a = l;ifIPLTPE = 0 

= 1.2; ifIPLTPE > 0 

b - l;ifISTATE-0 

2; if 1ST ATE > 0 

c = l;ifIPRINT = 0 

1.5; if IPRINT > 0 

d = l;ifIFDBAK = 0 

1.2; iflFDBAK > 0 
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^^FINDRMS "" 


(b»c*d) * 


NFREQ^NDOF^ 


(a^b^c^d) 

3000 


1000 J 
NSETS 

* £ NFREQ»(NLD|*NDOF) 

i= ! 


C^SORTOLS 

^^PpLOTQLS 

CPl221 


I second per item sorted 
1 second per item plotted 

10. + CPpijviDi^j^S + CPgQR-j-QLs + CPpLOTQLS 


Printed Output 

The printed output line limit has been set to 100000 lines which should be enough for 
any L221 (TEV156) execution. 


Punched Output 


When punched output is requested <1PUNCH = 1 on card 4.3 of FINDRMS), the number 
of cards punched will be 


T NSETS 

I NLD: 
' 1 


where: 


NSETS = 

Number of load sets 

NLDj 

Number of loads in set i 


6.3 CARD INPUT DATA 

The card input described in sections 6.3.2 to 6.3.5 is summarized in section 6.3.6. The 
summary is a quick reference for the necessary card input and is included for use only 
after familiarity with the program has been obtained. 

The taskis) performed by L221 (TEV156) are broken into three subtasks, each with its 
own section of code known as a primary overlay- The entire set of primary overlays is 
driven by a small program (main overlay) named L221. 

L221 reads program directive cards to: 

• Assure that the data being read is intended for L22I (TBV156) 
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9 Determine which section of code (primary overlay) of L221 ia to be executed next 

Each primary overlay has its own card input data sets. The primary overlays and data 
set numbers are listed in table 2. 


Table 2.— Primary Overlays and Data Set Numbers 


Overiay 

Data 

Purpose of primary overlay 

Number 

Name 

sets 

1 

FINORMS 

2 through 
49 

Calculate the generalized coordinates and 
loads over a range of frequencies and find 
the RMS values of loads. A, and the- 
number of zero crossings, N^. 

2 

SORTQLS 

[ 50 through 
61 

1 

1 Sort data generated by FI NORMS and write 
1 a magnetic file containing the data in 
vectors (frequency dependent) suitable for 
plotting. 

3 

PLOTQLS 

62 through 
79 

1 

Sort data generated by FIND RMS and write 
a magnetic file contair^rng plot instructions 
and data for the COMpSO plotter. 


The order in which 
displayed in figure 2. 


these card sets are input (and the overlays are executed) is 
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uJi] 



Figure 2.— Flow of L221 (TEV156) Card Input Data 
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CS.1 FORMAT OF CARD INPUT DATA 

card data are read in fixed fields-specific columns of the carda. The reqi^ived card 
columns are defined next to each keyword or variable on the pages that follow. Note the‘ 
following conventions used throughout the program: 

e All floating point variables are read with format ElO.O 

f All integer variables are read with format 15. 

^ All hollerith variables (keywords, etc.) are read with format AlO. 

Therefore, all data fields end on a card column which is a multiple of five. 

When the program is trying to recognize keywords, it checks only the first five 
characters. Any additional characters are ignored. 



Format of Matrices Input on Cards 

Matrices are read from cards by rows (each row begins on a new card) with the 
format 7E10.0. When a matrix contains only one column, the transpose of the matrix is 
read from cards. 

Example: 

If {Ml ] is to be input as a 9 by 9 matrix, 18 cards will be needed. 


L 

1 

Hid. 2) 1 

> • • 

1 

Ml (1.7) 1 

Card 1 1 

1 

TO 

• - 20 


60 

70 


L 

Ml0s8) 1 

Mid. 9) 1 


1 

_ 1 

Card 2 | 

1 

10 

20 



70 


L 

M,(2,1) 1 

Ml (2.2) 1 

■ • ■ 

1 

M,(2.7) I 

Card 3 | 

1 

10 

20 


60 

70 


L 

M,(2.8) 1 

Mi(2.9) 1 

■ • • 

1 

1 

Card 4 | 

1 

10 

20 

• 


70 





• 







• 




L 

Mt(9.1) I 

Mi(9.2) I 

• • ■ 

1 

Ml (9.7) 1 

Card 17 | 

1 

10 

20 


60 

70 


L 

M,(9.8) 1 

Ml (9.9) 1 

• • • 

1 

1 

Card 18 | 

1 

10 

20 



70 



31 



6.3.2 PROGRAM L221 


The driving module of L221 (TEV156) (the main overlay is known as L221) will read 
and check the first five characters of the input card to determine the type of execuition 
that is to take place. 

Note: All underlined capital characters contained in the KEYWORDA^ARIABLE field 
of the input card sets must be left justified and punched in the card columns 
specified in the COLS, field of the input card sets. 

Card Set 1.0-Introduce L221 (TEV156) Card Input Data 


Card 1.1-Frequency Response Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

$FREQuency 

AlO 

This card must be the first card read by L221 (TEV156). It signifies 
that the input card data is to be used by L221 (TEV156) for a Fre- 
quency Response Analysis. 

11-70 



Available for comments. 


Card 1.2-Title Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

$TITLe 

AlO 

A title card that Is Intended to label the program's printed output 
with miscellaneous titling Information. 

After this card is read and printed. control is returned to the 
main overlay for the reading of the next instructional keyword card. 
The number of $TITLe cards Is unlimited. 

11-70 



Available for comments. 


Card 1.3-Core Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

$C0RE 

AlO 

A keyword designating a card which sets the maximum field length to 
be used by L221 (TEV156). 

11-20 

MAXCOR 

010 

The maximum field length, in octal, which L221 (TEV156) may use. 
Default: 200000g 

21-70 



Available for conments. 


32 




Card 1.4-Special Programmer Checkout Option 


COLS. 

JCEYMORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

i-id“ 

[ICHECk 

Aid' 

This card is used as an aid in the checkout of fatal error diag- 
nostics. This option will allow L221 (TEV156) to continue execu- 
tion even after a fatal error has been diagnosed. 

11-70 



Available for comments. 


To use this option, the programmer must use the following card sets: 

Card 1.1 $FREQuency 

Card 1.4 $CHECk 


Any card sets the programmer wishes to be tested. 


Card set 80.0 


$QUIT 




641^ PROGRAM FINDRMS 


The card input data for FINDRMS are divided into three main blocks: 


BLOCK I 


Options and constant data 


BLOCK II Instructions and matrices needed to calculate the generalized 

coordinates frequency response 

BLOCK III Instructions and matrices needed to calculate the loads frequency 

response 


The three blocks are always in this order, although BLOCK III may be repeated as 
many as 10 times. BLOCK I is to be input on cards, but the main parts of BLOCKS II 
andlll-the matrices-may be input on cards or on the magnetic files "INTAPE” 
and "LDTAPE”. 


The number of 
sets (NSETS) 
is limited to 10. 


( Entry intoA 
FINDRMS J 


BLOCK I 



^ Return to A 
V ^Main Overlap 


Card Sets 2-11 


Card Sets 12-30 


Card Sets 31-48 


Card Set 49 


BLOCK I-Options and Constant Data 

The $FINDRMS card (card set 2.0) introduces BLOCK I, card input data defining the 
necessary options and constants. The flow of BLOCK I data cards is displayed in 
figure 3. 





Figure 3.— Flow of FINDRMS Block / Card Input Data 
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Card Set 2.0-Call Card for FINDRMS 


1 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SFINDrms 

AlO 

This Card Set Initializes the call to the primary overlay* FINDRMS. 
FINDRMS win read Card Sets 3-49 and perform the necessary calcu- 
lations. Control win be returned to the main overlay* L221 upon 
the reading of Card Set 49. 

^11-70 



Available for comments. 


Card Set 3.0-Label Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-70 

LABEL 

7A10 

A one card data case label which Is printed above the resulting 
output. 
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Cilrd ^et and File Names 

Card 4.1-Probliem Siie 


COLS. 

. . j 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1^5 

NDOF 

IS 

The number of degrees of freedom, I.e. the number of generalized 
coordinates or equations of motion. 

(0 < NDOF < 100) 

— 

6^10 

NFREQ 

IS 

The number of frequencies over which the solution will be found. 
(HFREQ >2) 

11-15 

NPAN 

15 

The number of panels for gradual penetration. If the gradual 
penetration option Is not desired NPAN must be set to zero. 

When using DYLOFLEX tapes for “ INTAPE" and/or "LDTAPE" the nuab4r 
panels must be >1 because gradual penetration Is always assiMSd 
In DYLOFLX, thus no {C 3 }. 




(0 < NPAN S 100) 

16-20 

NKVAL 

15 

The number of k-values used in an unsteady aerodynamic analysis. 

If NKVAL is greater than zero, matrices will be read at the defined 
frequencies (Card Set 17) and Interpolation will occur to find the 
matrices at the solution frequencies. IC 3 ) or [^] are complex with 
NKVAL > 0. When NKVAL - 0 all input matrices are real. When 
using DYLOFLEX tapes for "INTAPE" and/or "LDTAPE" the number 
of k values must be > 1 . 

(NKVAL > 0) 

Default: 0 




21-25 

NFORC 

15 

The number of excitation forces for the handling of multiple 
forcing functions (see Appendix C). 

Note; If more than one excitation force Is being used the fol- 
lowing limitations are Imposed: 

• NPAN = 0. 

• no feedback. 

• no loads correlation. 

• the frequency responses and load transfer functions may 
be plotted for only the first excitation forces. 

• IC 2 } and j € 2 } are null. 

• the loads per degree of freedom can not be printed, 
(NFORC > 1) 

Default: 1 

26-30 

NDOFD 

15 

The number of degrees of freedom to be deleted. 

Note: See Card 7.2 for the input of the degrees of freedom to be 

deleted. 

(NDOFD < NOOF) 




Card 4.2-Optional Capabilities 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

I FREQ 

15 

This option chooses the type of frequency matrix Input. 

« -1 The frequencies will be read from cards In cycles per second, 
(see Card 11.1). 

* 0 The frequencies will be read from cards In radians per second, 
(see Card 11.2). 

> 1,2, 3,4, 5 The frequency matrix will be generated according to 
Appendix A, using an Initial frequency, FREQ^» and an in- 
crement DELTAF In cycles per second, (see Card 11.3). 

6-10 

( 

ISPEC 

15 

This option specifies the type of Input spectrum to be used, 

Cn Is in radians per unit distance). 

» 0 Oryden . o 

[1 + (Ln)^]^ 

= 1 Von Kaniian - 

L 1 + 8/3 (1.339Ln)^ 

♦(n) - 5 -rr>^ 

[1 + (1.339L£ir]'''® 

= 2 Exponential 

4>{n) =• K (n)*” 

Note: See card Set 6 for L and Card 8.1 for ic and n. 

= -NSP 

A tabular Input spectrum defined at NSP points (see Cards 8.2, 
8.3, and 3.4). 

ii-15 

_ — ^ — 

I STATE 

15 



This option requests a static-elastic solution. 

- 0 No static-elastic solution. 

« 1 A static-elastic solution will be found after zeroing columns 
3 through NDOF of [Mg], [Mj], [M 5 ]. [Hg], [Sj], [B 3 ]. [B 5 ]. 
and [Hg]. Implying that the rigid body degrees of freedom are 
In columns 1 and 2 . 

« 2 A static-elastic solution will be found, by zeroing columns of 
matrices [M 2 ]. [Mj], [Mg], [Mg], [Rg]. [B 3 ]. [Rg], and [1^] 
according to Card Set 7. 
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16-20 

IFDBAK 

15 

This option requests that a feedback loop be used. 

* 0 No feedback loop used, 

= 1 A standard feedback loop will be used. 

* 2 A special feedback loop will be used, (see RLIMIT of Card 

10.1 and Cards 10.6, 10.7, and 10.8). Also, see Appendix B. 

21-25 

IDAMP 

15 

This option specifies the number of different structural damping 
factors. 

= 0 The same factor will be used for all freedoms, (see Card 9.1). 
= 1 A different factor will be used for each frequency, (see Card 
9.2). 

26-30 

INTZRO 

15 

This option chooses the starting point for the RMS Integration. 

* 0 The RMS integration will start at zero. The output spectrum 
is assumed to grow linearly between the frequencies 0 and 
FREQj (from the value 0. to the value calculated at FREQ^). 
The area in the resulting triangle is Included in the 
integration. 

- 1 The RMS Integration will start at FREQj. The range between 
frequencies 0 and FREQj is Ignored. 
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I 


Card 4.3-Output Options 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

IPRINT 

15 

This option chooses the type of printing to take place. 

* 0 Print only the final T values and the nunber of zero 
crossings. 

- 1 Print the generalized coordinates and load transfer functions 
for each frequency as well as the final T and number of zero 
crossings. 

6-10 

IRMSPR 

15 

This option chooses the type of printing for theT values. 

= 0 Print only the final A values. 

> 0 Print the A values for each of the last IRMSPR frequencies. 
(IRMSPR < NFREQ) 

11-15 

I PUNCH 

15 

This option requests that the final A values and nunber of zero 
crossings be punched on cards. 

* 0 No punched output. 

= 1 Punch on cards the final A values and nunber of zero 

crossings. 

16-20 

ICKPRT 

15 

This option chooses the amount of checkout printout for the problem. 
= 0 The input matrices are not printed. 

= 1 The input matrices will be printed as read. 

= 2 The input matrices plus intermediate data for the first 

three frequencies will be printed. 

= 5 Special program checkout feature. It prints 1 and 2 above 

and a record of all VARDIM routines called. 

21-25 

IPLRMS 

15 

This option requests the writing of final "A values on the plot 
data file “IPLTPE" 

= 0 7 values n^ saved, 

• 1 T values will be written on "IPLTPE" for each frequency, 
(see Card 4.4). 
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26-30 

LPDOF 

15 

This option requests the printing of loads per degree of freedom 
at certain frequencies. 

= 0 The loads per degree of freedom will not be found. 

> 0 The loads due to each degree of freedom will be found and 




printed for the LPDOF frequencies specified on Card 11.4. 




(LPDOF 5 10) 


Card 4.4-File Names 

Note: All magnetic files not to be defined must have their appropriate fields left 

blank. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

INTAPE 

AlO 

This option indicates the location of the input matrices. 

- The input matrices will be read from cards if the field is left 
blank. 

- The input matrices will be read from a magnetic file in READTP/ 
WRTETP format. The file name must be a left-justified 1 to 7 
character name with the first character being alphabetic. 

11-20 

IPLTPE 

AlO 

This option indicates whether or not a plotting magnetic file will 
be created. 

- The plotting information will not be written if the field is 
left blank. 

- The plotting information will be written on a magnetic file in 
REATP/WRTETP format. The magnetic file name must be a left- 
justified, 1 to 7 character name with the first character being 
alphabetic. 

21-30 

IRTAPE 

AlO 

This option indicates whether or not the generalized coordinates and 
the load transfer functions will be written on a magnetic file. 

- The generalized coordinates, and the load transfer functions will 
not be written onto a magnetic file if left blank. 

- The generalized coordinates, and the load transfer functions will 
be written on a magnetic file in READTP/WRTETP format. The 
magnetic file name must be a left-justified, 1 to 7 character 
name with the first character being alphabetic. 
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31-40 

IFTAPE 

AlO 

This option Indicates the location of the Tabular Input Spectra. 

- The Tabular Input Spectra will be read from cards If the field 
Is left blank. 

- The Tabular Input Spectra will be read from a magnetic file In 
READTP/WRTETP format. The file name must be a left justified 

1 to 7 character file name with the first character being 




alphabetic. 


Card Set 5.0-Wagner and Kiassner Indicial Lift Growth Function Coefficients 
Card 5.1- Wagner Indicial Lift Growth Coefficients 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

^1 

ElO.O 

The Wagner Function Coefficients. 

11-20 


ElO.O 

Bi|S C-iS D<iS 

g(s) * s+B^ " s+Y] 

21-30 


ElO.O 

31-40 


ElO.O 

where s =1(u and w Is the frequency In radians per 
second. 

(See table 3) 

41-50 

“1 

ElO.O 

51-60 


ElO.O 

Note; When NKVAL > 1 of card 4.1, set Ai - 1. and 
all other variables to zero. 

61-70 


ElO.O 


Card 5.2-Kussner Indicial Lift Growth Function Coefficients 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

*2 

ElO.O 

The Kussner Function Coefficients. 

11-20 

^2 

ElO.O 

8,s C,s D,s 

f(s) . Ao- . _2_ . 2 

21-30 

^2 

ElO.O 

^ ' 2 s+og S+Cg S + Y 2 

31-40 

'>2 

ElO.O 

Where s =103 and w is the frequency in radians per 
second. 

(See table 4) 

41-50 

«2 

ElO.O 

51-60 

h 

ElO.O 

Note: When NFORC > 1 or NKVAL > 1 set A 2 - 1. and 
all other variables to zero. 

61-70 

^2 

ElO.O 


Tables 3 and 4 list the different values for the Wagner and Kussner functions as 
functions of aspect ratio, AR, at M = 0. These values were obtained from reference 3. 
The user is urged to refer to reference 3 and 4 in order to obtain a greater 
understanding of the derivation of these indicial functions and the theoretical aspects 
associated with their derivation. 

Table 3.— Coefficients for the Wagner indicial Lift Growth Function at M = 0 


AR 



Cl 

Dl 

“1 

h 

1-1 

QO 

1.0 

.165 

.335 

0 

.045 

.30 

0 

6 

1.0 

.361 

0 

0 

.381 

0 

JL. 

3 

1.0 

.283 

0 

0 

.540 

0 

0 


Table 4.— Coefficients for the Kussner Indicial Lift Growth Function atM = 0 
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Card Set 6.0-Input Constant 


COLS. 

1-io 


ii -20 


KEYWORD/ 

VARIABLE 

Z 


VEL 


FORMAT 

ElO.O 


ElO.O 


DESCRIPTION 

The scalar constant for the equations of motion and loads equa- 
tion convoluted forcing function coefficient matrices, (e.g. Ic^} 
and {C 3 }). In DYLOFLX, Z must be input as a negative number to 
obtain sign compatibility for the gust forcing function generated 
in L217(E0H) and L213(LDADS) 

(Z/0) 

The velocity (true airspeed) of the frequencies are converted to 


radians per unit distance by the following conversion: 


where n is in radians/distance and w is in radians per second. 


21-30 


L 

(TSCALE) 


ElO.O 


is used to determine the input spectrum.* 

(VEL > 0) 

The turbulence scalar constant. It is required when ISPEC of 
Card 4.2 is equal to 0 or 1 . (See eqn 26 and eqn 27 Of Section 
5.0 for the Dryden and Von Karman input spectra). 


31-40 


C 

(GPSCAL) 


ElO.O 


Note: The unit distance must be in the same units described in 

VEL above. 

(L > 0) 

The gradual penetration scalar constant. (See eqn 19 of Section 5.0 
for the calculation of gradual penetration). 

Note: This variable is needed only when gradual penetration is 

requested. 



ElO.O 


(C^O) __ 

The scaling factor for the frequencies at which the Tabular Input 
Spectra are defined. (See Card 8.2). 

(Default = 1) 


•For DYLOFLEX when gradual penetration is used, the length unit of VEL must be the same as 
the f units of length established in L217 (EOM). 
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51-60 

TABS 

£10.0 

The scaling factor for the Tabular Input Spectra. ^ee Card 

8.3 or 8.4). 

(Default ■ 1) 

61-70 

FDSCAL 

£10. 0 

The scalar constant for the frequencies contained In |FR£QM}, 
the array of frequencies at which Input matrices are defined. 

This scalar can be used only when NKVAL of Card 4.1 is greater 
than zero. 

Note: If the frequencies are to be read in as radians /second, 

and {FR£QM| Is In cycles/second, then FDSCAL must = 
6.2831853, (which Is 27r). 

If the frequencies are to be read In as cycles/second, and 
{fR£QM} is In radians/second, then FDSCAL must » 
1/6.2831853, (which Is 1/2tt). 

(Default = 1) 
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Card Set 7.0-Deletion of Degrees of Freedom 
Card 7.1-Deletion of Freedoms for a Static-Elastic Solution 
Omit card 7.1 if ISTATE ^ 1 on card 4.2. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

NSTELj 

1415 

NSTEL^-, 1=1, NDOF where NDOF (Card 4.1) is the number 
of degrees of freedom, 

NSTEL chooses the columns to be zeroed in the matrices 

6-10 

NSTEL 2 


[M 2 ]. [M 3 ]. [Mg], [Mg], [M 2 ]. [H 3 ], [Mg], and [Hgjwhen performing 
a static-elastic solution. 

• 

• 

• 


• 

• 


If NSTEL i = 0, then column 1 of the above matrices will be zeroed. 


NSTELnoof 


Repeat the card with 14 numbers per card until all elements are 




defined. 


Card 7.2-Deletion of Degrees of Freedom 
Omit card 7.2 if NDOFD = 0 on card 4.1. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

NOEL^ 

1415 

NDEL-j, i = 1, NDOFD where NDOFD (Card 4.1) is the number of degrees 

6-10 

NDEL 2 


of freedom to be deleted. 

11-15 

• 

• 

NDEL 3 

• 

• 


NDELi is the specified degree of freedom to be deleted. 

(NDELi <ND0F) 

For the generalized coordinates, the degrees of freedom will be 
deleted by zeroing the appropriate rows and columns of [SUMM] 
and the rows of |SUMC| . Also, the diagonal elements of the 

• 

• 


appropriate rows of [SUMM] will be set equal to 1. 




In the loads solution the degrees of freedom will be deleted by 
zeroing the appropriate columns of [SUMM] . 


”°^*-ND 0 F 0 


Repeat the card with 14 numbers per card until all elements are 
defined. 




Card Set 8.0-Input Power Spectra 


Omit card set 8.0 if ISPEC = 0 or 1 on card 4.2. 
Card 8.1 “Exponential Power Spectrum Input 
Read card 8.1 if ISPEC = 2 on card 4.2. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

K 

(OCOEF) 

ElO.O 

The input power spectrum coefficients. 
i)>(n) = 

(See eqn 28 of Section 5.0 for the exponential input spectrum). 

11-20 

V 

(OEX) 

ElO.O 


Card 8.2-Frequencies of the Tabular Input Spectrum 

Omit card 8.2 if ISPEC ^ 0 on card 4.2 or if IFTAPE of card 4.4 is defined. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

FREQS 1 

7E10.0 

FREQS^ , i = 1, NSP 

11-20 

• 

• 

• 

^ ... 

cr 

LiJ 

a: 

u. 


where NSP of Card 4.2 is the value for the Tabular 
Input Spectrum. 

The Tabular Input Spectrum Frequencies in cycles per seconds. 


FREQS, sp 


Repeat the card with 7 numbers per card until all elements are 
defined. 


FREQSi are the frequencies at which the tabular input spectra are defined. The 
frequencies must cover the range of solution frequencies defined in card set 11.0 and be 
defined in the same units. The scalar TABFS of card set 6.0 may be used to provide 
compatibility. L221 (TEV156) will not extrapolate. 
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Card 8.3-Single Tabular Input Spectrum 

Omit card 8.3 if ISPEC > 0 on card 4.2, or if NFORC > 1 on card 4.1, or if IFTAPE 
of card 4.4 is defined. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

PHIOSj 

7E10.0 

PHIOSi ; 1 » 1, NSP 

The value of the tabular input spectrum at FREQS{ ; i ° 1, NSP 
where NSP » (lSPEC|from Card 4.2. 

11-20 

PHI0$2 

• 

• 

• 

• 


PHI0S,5p 


Repeat the card with 7 numbers per card until all elements are 
defined. 


PHIOSi is a real number (no imaginary part). The scalar TABS in card set 6.0 may be 
used to modify the units. The program will interpolate to find the spectrum value at 
solution frequencies. 


1 
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Card 8.4-Multiple Tabular Input Spectrum 

Omit card 8.4 if ISPEC > 0 on card 4.2, or if NFORC 1 on card 4.1, or if IFTAPE of 
card 4.4 is deAned. 

Card 8.4 must be repeated for all combinations of (I,J) where 1=1, NFORC and J = I, 
NFORC. Because of symmetry, only the upper right triangulsu* combinations are input. 
Input row 1 first (J = 1, NFORC), then row 2 (J = 2, NFORC), . . ., and, finally, row 
NFORC (J = NFORC). 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

PHI0S(I,J)j 
(real ) 

7E10.0 

PHIOS(I.J),^; k = 1,NSP 

The value of the (1,0) tabular input spectrum at FREQS ; 
k =1, NSP where NSP = |lSPEC) from Card 4.2. 

Note that each element is a complex number with real and 
imaginary parts. 

11-20 

PHIOS(I.J)^ 

(imaginary) 

21-30 

PHI0S(I,J)2 
(real ) 

31-40 

PHI0S(I,J)2 

(imaginary) 

• 

• 


PHI0S(I,J)f,5p 
(real ) 


Repeat the card with 7 numbers per card until all elements are 
defined. 


PHI0S(I,J)^5P 

(imaginary) 




Card Set 9.0-Structural Damping 
Card 9.1-Same Damping Factor for All Freedoms 
Read card 9.1 if IDAMP = 0 on card 4.2. 



COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 


1-10 

9sd, 

ElO.O 

A structural damping factor used for all degrees of freedom. 


Card 9.2-Different Damping Factors for Each Freedom 
Read card 9.2 if IDAMP = 1 on card 4.2. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

5'sD, 

7E10.0 

9sd.. i * 1. NDOF 

where NDOF of Card 4.2 is the number of degrees of freedom. 
The structural damping factor for each degree of freedom. 

11-20 


• 

• 

• 

• 

• 


Hoof 


Repeat the card with 7 numbers per card until all elements are 
defined. 


63 




Card Set 10.0-Feedback Coefficients 


Omit card set 10.0 if IFDBAK = 0 on card 4.2. 
Card 10.1-Standard Feedback Coefficients 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

G 

ElO.O 

Coefficient of the feedback transfer function, FYDAMP, (see the 
equation below). 

11-20 

RLIMIT 

ElO.O 

The rudder limiting angle in radians. This variable is needed only 
when IFDBAK of Card 4.2 is equal to 2. 


Card 10.2-Standard Feedback Coefficients 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

Pi 

7E10.0 


11-20 

P2 


1 =1,7 

• 

• 

• 

• 


Coefficient of the feedback transfer function, FYDAMP, (see the 

• 

• 



61-70 



equation below). 


Card 10.3-Standard Feedback Coefficients 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 


7E10.0 


11-20 



U^, i = 1 .7 

• 

• 


Coefficient of the feedback transfer function, FYDAMP, (see the 

• 

% 



• 

• 


equation below). 

61-70 

^7 
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Card 10.4-Standard Feedback Coefndents 


■ 

COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

MO 


7E10.0 

R^. 1 ■ 1,7 

Coefficients of the feedback transfer function, FYDAMP, (see the 
equation below). 

11-20 


• 

• 

m 

• 

m 

• 

61-70 

«7 


Card 10.5-Standard Feedback Coefficients 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

MO 


7E10.0 

T^, i » 1,7 

Coefficients of the feedback transfer function, FYDAMP, (see the 
equation below). 

11-20 

^2 

* 

• 

• 

• 

^1-70 

T 7 


> Fydamp = 


G(P|S^ + Uis + Rj)(P28^ + U28 + R2)(P38^ + U38 + R3KT18 + 1 ) 
(P38^ + U38 + R3)...(P78^ + U78 + R7KT2S + IKT38 + 1)...(T78 + 1) 
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Omit cards 10.6 through 10.8 if the standard feedback loop is chosen (IFDBAK « 1 on 
card 4.2). 

Card 10.6-Special Feedback Coefficients (see app. B for coefficient definitions) 


COLS. 

keywo.ro/ 

VARIABLE 

FORfW 

DESCRIPTION 

1-10 

CS2j 

7E10.0 

CS2^, i * 1. NDOF 

where NDOF (Card 4.1) is the number of degrees of freedom. 

11-20 

CS22 

• 

* 

• 


^^^NDOF 


Repeat the card with 7 numbers per card until all elements are 
defined. 


Card 10.7-Special Feedback Coefficients (see app. B for coefficient definitions) 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

CSl^ 

7E10.0 

CSl^ 1 = 1, NDOF 

where NDOF (Card 4.1) is the number of degrees of freedom. 

11-20 

CS 12 

J 

• 


*'^^N00F 


Repeat the card with 7 numbers per card until all elements are 
defined. 


Card 10.8-Special Feedback Coefficients (see app. B for coefficient definitions) 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

CSOj 

7E10.0 

CSO^, i = 1, NDOF 

where NDOF (Card 4.1) is the number of degrees of freedom. 

11-20 

CS 02 

• 

• 

• 



‘•^^NDOF 


Repeat the card with 7 numbers per card until all elements are 
defined. 
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Card Set 11.0-Frequency Array 

Card 11.1 -Frequencies in Cycles per Second 
Read card 11.1 if IFREQ =-l on card 4.2. 


COLS. 

KEYWORD/ 

VARIABLE 

1-10 

' freqY 

11-20 

FREQg 

• 

• . 

'"’^^‘^NFREQ 


FORMAT 

7E10.0 


DESCRIPTION 

FREQ^. , i = 1, NFREQ 

where NFREQ (Card 4.1) is the number of frequencies. 

The input frequencies are in cycles per second. 

Repeat the card with 7 numbers per card until all elements are 
defined. 


Card 11.2-Frequencies in Radians per Second 
Read card 11.2 if IFREQ = 0 on card 4.2. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

FREQj 

7E10.0 

FREQ., 1=1, NFREQ 

where NFREQ (Card 4.1) is the number of frequencies. 
The input frequencies are in radians per second. 

Note: Internally the input frequencies are converted 

from radians per second into cycles per second. 

11-20 

FREQg 


'"‘^^‘^NFREQ 

Repeat the card with 7 numbers per card until all elements are 

defined. 

. .. 
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Card 11. 3-First Frequency and an Increment 
Read card 11.3 if IFREQ > 0 on card 4.2. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

FREQj 

ElO.O 

The first frequency (cps). 

11-20 

OELTAF 

ElO.O 

The frequency increment in cycles/ second. 


Card 11.4-Frequencies at Which the Loads per Degree of Freedom Will Be Printed 

Read card 11.4 if LPDOF > 0 on card 4.3. 

Note: If IFREQ = 0 on card 4.2, FREQLi must be in radians per second. 

If IFREQ = -1 on card 4.2, FREQLi must be in cycles per second. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

FREQL^ 

7E10.0 

FREQL^. i = 1, LPDOF 

11-20 

FREQL 2 


where LPDOF is from Card 4.3. 

• 

’ 

* 


The frequencies at which the loads per degree of freedom 


• 


will be printed. 


FREQLlpjjjjp 


Repeat the card with 7 numbers per card until all elements 




1 are defined. 


Note: Given FREQLj, FINDRMS will choose the element of the array {FREQ} which is 
equal to or greater than FREQLj. 

FREQL| must be ^ FREQnfrEQ* 
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BLOCK Il-Instructioiia and Matrices Needed To Find the Generalised 
Coordinates 

BLOCK II, card sets 12.0 through 30.0, always follows BLOCK I. The flow of BLOCK II 
data cards is displayed in figure 4. BLOCK II is always introduced by card set 12.0. 

Card Set 12.0-Introduce the Equations of Motion 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

$EqUAt1on 

AlO 

This card introduces the cards defining the equations of motion. 

11-70 



Available for coinnents. 





Figure 4.— Flow of F I NORMS Block II Card Input Data 










2A 



* Indicates that card 
is optional depending 
upon variables on 
Cards 13.2 and 13.3. 


NK = 1 .NKVAL 



Figure 4.— (Continued) 


Loop on number 
of frequencies 
* at which Input 
matrices are 
defined. 
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[$3 

Card Set 30* 


<NKVAL. 


= NKVAL 


* Indicates that card 
is optional depending 
upon variables on 
Cards 13.2 and 13.3 


Figure 4.— (Concluded) 



I 


Card Set 13.0-Choose the Non*Null Matrices 

Card set 12.0 must be followed by either card 13.1 or 13.2. Card 13.1 indicates that the 
equations of motion matrices are on magnetic file "INTAPE” which begins with a 
DYLOFLEX header array (see sec. 6.4). Card 13.2 must be used if the matrices are on 
cards or if "INTAPE” does not contain the DYLOFLEX header array. 

Card 13.1-Signal an "INTAPE” With DYLOFLEX Header Array 

Omit card 13.1 if the equations of motion are on cards or if "INTAPE” does not 
contain the DYLOFLEX header array. 


COIS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

MO 

OVLOFLX 

AlO 

Keyword Indicates that ** INTAPE ** begins with a DYLOFLEX header array. 

11-70 



Available for comments. 


|NTAPE (card 4.4) must be defined in order to use card 13.1. When card 13.1 is used, all 
other BLOCK II cards (13.1 through 30.0) must be omitted. The DYLOFLEX header 
array and the equations of motion will be read from magnetic file "INTAPE”. The 
header array at the beginning of "INTAPE” specifies the matrix sizes and indicates 
which ones are not null and appear on the file (see sec. 6.4) 
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Card 13.2-Choose Non-Null Equations of Motion Matrices 
Omit cards 13.2 and 13.3 if card 13.1 was used. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

NMl 

15 

NMl. 1 » 1,6 

6-10 

NM2 

15 

If NMl = 0 the corresponding equation of motion matrix [M^] 

11-15 

NM3 

15 

will not be read. It is assumed to be null. 

16-20 

NM4 

15 


21-25 

NM5 

15 

If NMi > 0 the matrix [Mi] will be read. 

26-30 

NM6 

15 


31-35 

NC2 

15 

If NC2 = 0 the corresponding equation of motion forcing function 
coefficient matrix jC 2 ^ or [Cg] will not be read. 
It is assumed to be null. 

If NC2 > 0 the matrix {C 2 } or [C 2 ] will be read. 

36-40 

NC3 

15 

If NC3 = 0 the corresponding equation of motion convoluted forcing 
function coefficiet matrix |C^| or [C^] will not be 
read. It is assumed to be null. 

If NC3 > 0 the matrix jC^I or [C^] will be read. 

Note: If the number of panels (NPAN of Card 4.1} is > 0. then 

NC3 is set to zero. 

41-45 

NFLL 

1 

15 

If NFLL = 0 the corresponding matrix of the streamline distance In 
feet from the point first encountering gusts to the 
points encountering gusts later will not be read. It 
is assumed to be null. 

If NFLL = 1 the matrix If^} will be read. 

Note: If the number of panels (NPAN of Card 4.1) is « 0 then 

NFLL should be zero.. 


64 




46-50 

NPHIT 

15 

If NPHIT = 0 

the corresponding matrix of the gust forces acting on 





each gradual penetration lifting panel will not be read. 




If NPHIT = 1 

the matrix [$] will be read. 




Note: If the number of panels (NPAN of Card 4.1) is = 0, then 




NPHIT 

should be zero. 


Card 13.3-Choose Non-Null Feedback Matrices 

Omit card 13.3 if card 13.1 was input or if IFDBAK of card 4.2 = 0. 


COLS. 

1-5 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

NSl 

15 

NSi. 1-1,6 

= 0 the corresponding coefficient of feedback for the equations 
of motion matrix, [S^] , will not be read. It will be assumed 

null . 

> 0 the corresponding matrix, [Si] , will be read. 

6-10 

NS2 

15 

11-15 

NS3 

15 

16-20 

NS4 

15 

21-25 

NS5 

15 

26-30 

NS6 

15 


Omit card sets 14.0 through 30.0 if INTAPE of card 4.4 is deflned. Instead, the matrices 
will be read from the magnetic file "INTAPE”. 

Card Set 14.0 [M^] -Generalized Stiffness Matrix 

Omit card set 14.0 if NMl = 0 on card 13.2. [Mi ] is of size NDOF x NDOF. 

Card Set 15.0 [M 2 ] -Viscous Structural Damping Matrix 

Omit card set 15.0 if NM 2 = 0 on card 13.2. [M 2 ] is of size NDOF x NDOF. 

Card Set 16.0 [M 3 ] -Generalized Mass Matrix 

Omit card set 16.0 if NM3 = 0 on card 13.2. [M 3 ] is of size NDOF x NDOF. 

Card Set 17.0 (FREQM [-Frequencies at Which the Aerodynamic Equations of 

Motion Are Defined 

Omit card set 17.0 if NKVAL = 0 on card 4.1. {FREQM} is of size 1 x NKVAL. 
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Repeat card sets 18.0 through 30.0 NKVAL TIMES: NK = 1 , NKVAL. 

Card Set 18.0 [M 4 ] -Convoluted Generalized Coordinate Displacement Matrix 

Omit card set 18.0 if NM4 = 0 on card 13.2. [M 4 ] is of size NDOF x NDOF. 

Card Set 19.0 [M 5 ] -Convoluted Generalized Coordinate Rate Matrix 
Omit card set 19.0 if NM5 = 0 on card 13.2. [M 5 ] is of size NDOF x NDOF. 

Card Set 20.0 [M 0 ] -Convoluted Generalized Coordinate Acceleration Matrix 
Omit card set 20.0 if NM 6 = 0 on card 13.2. [Mg] is of size NDOF x NDOF. 

Omit card sets 21.0 through 26.0 if IFDBAK = 0 on card 4.2. 

Card Set 21.0 [Si]-Feedback Coefficient Matrix for Generalized Coordinate 

Displacement 

Omit card set 21.0 if NSl = 0 on card 13.3. [Sj] is of size NDOF x NDOF. 

Card Set 22.0 [S 2 ] -Feedback Coefficient Matrix for Generalized Coordinate Rate 

Omit card set 22.0 if NS 2 = 0 on card 13.3. [S 2 ] is of size NDOF x NDOF. 

Card Set 23.0 [SsJ-Feedback Coefficient Matrix for Generalized Coordinate 

Acceleration 

Omit card set 23.0 if NS3 = 0 on card 13.3. [S 3 ] is of size NDOF x NDOF. 

Card Set 24.0 [S 4 ] -Feedback Coefficient Matrix for Convoluted Generalized 

Coordinate Displacement 

Omit card set 24.0 if NS4 = 0 on card 13.3. [S 4 ] is of size NDOF x NDOF. 

Card Set 25.0 [Ssj-Feedback Coefficient Matrix for Convoluted Generalized 

Coordinate Rate 

Omit card set 25.0 if NS5 = 0 on card 13.3. [S 5 ] is of size NDOF x NDOF. 

Card Set 26.0 [SeJ-Feedback Coefficient Matrix for Convoluted Generalized 

Coordinate Acceleration 

Omit card set 26.0 if NS 6 = 0 on card 13.3. [Sg] is of size NDOF x NDOF. 

Card Set 27.0 {C 2 }-Forcing Function Coefficient Matrix 

Omit card set 27.0 if NC 2 = 0 on card 13.2. {C 2 } is of size NDOF x 1 . Input as a row. 
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Card Set 28.0 {C 3 } -Convoluted Forcing Function Coefficient Matrix 

Omit card set 28.0 if NC3 = 0 on card 13.2 or if NPAN > 0 on card 4.1. 

If there is a single forcing function (NFORC = 0 on card 4.1), {C3} is of size NDOFx 1 
and is input as a row. 

However, if the equations of motion are frequency dependent (NKVAL > 0 on card 4.1), 
{C3} is complex and must be input as a row of size 2*NDOF with alternating real and 
imaginary parts. 

If there are multiple forcing functions (NFORC > 0 on card 4. 1). [C3] is of size 
NDOF X NFORC. 

However, if in addition the equations of motion are frequency dependent (NKVAL > 0 
on card 4.1), [C3] is complex and must be input by rows with 2*NFORC numbers per 
row, alternating real and imaginary parts. 

Card Set 29.0 {f j^l-Distances Between First Panel and Other Gradual 

Penetration Panels 

Omit card set 29.0 if NPAN = 0 on card 4.1. 

Also, card set 29.0 is input only for the first frequency, {f^} is of size NPAN x 1. 
Input as a row. 

Card Set 30.0 [c/>]-Gust Forces Acting on Each Gradual Penetration Lifting 
Panel 

Omit card set 30.0 if NPAN = 0 on card 4.1. [^] is of size NDOF x NPAN. 

[<^] contains real numbers if NKVAL = 0 on card 4.1. 

[0] contains complex numbers if NKVAL > 0 indicating that the matrices are 
frequency dependent. When [0] is complex, it is input with alternating real and 
imaginary parts; 2* NPAN numbers to a row. 


BLOCK Ill-Instructioiis and Matrices Needed for the Loads Solution 

BLOCK III, card sets 31.0 through 48.0, follows BLOCK II. BLOCK III may be repeated 
10 times, once for each load set. The flow of BLOCK III data cards is displayed in 
figure 5. BLOCK III is always introduced by card set 31.0. The end-of-load sets and the 
end of FINDRMS execution are indicated by card set 48.0, $END. 

Card Set 31.0-Load Instructions 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

$L0AD 

AlO 

The keyword indicating that a new load set will be processed. 

11-15 

NLO 

15 

The number of load equations In this load set. 
(NLO < 100) 

16-20 

ICORR 

15 

The number of loads to be correlated with other loads. If NFORC 
of Card 4.1 >1, loads correlation can not be used. 

(ICORR < NLD) 

21-30 

c 

(GPSCLB) 

ElO.O 

The scaling factor which is only used for gradual penetration. 
See eqn 22 of Section 5.0. 

(C > 0) 

31-40 

LDTAPE 

AlO 

The magnetic file name that will contain the load equations input 
matrices. 

The tape name is defined as a one to seven character name, 
beginning with an alphabetic character and being left justified. 

If LDTAPE is left blank it will default to being the same as INTAPE 
on Card 4.4. 

If the load equations are to be input by cards, LDTAPE is left blank. 

Note: Load equation matrices can only be 
read in by cards if the equations 
of motion matrices are read in 
by Cards, that is, if INTAPE 
on Card 4.4 is left blank 













Loop on number 
of frequencies 
at which Input 
matrices are 
defined. 


Figure 5.— (Continued) 
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Card Set 32.0-Load Correlation Specifications 

r; ’Omit card set 32.0 if ICORR of card set 31.0 = 0, or if NFORC of card 4.1 > 1. . 
Repeat card set 32.0 ICORR times. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

LOAD^ 

15 

The load number that is to be correlated with up to 10 other 
loads. 

6-10 

LOADC(^^l) 

15 

The loads to be correlated with LOAD^ . 
(LOAD^ < NLD) 

11-15 

L0ADC(i^2) 

• 

• 

• 

• 

• 

• 

51-55 

LOADC(i^lO) 


The load number LOADj will be correlated with the loads numbered (LOADCqj), 
j.= 1, 10). 



Card Set ^3.0*>Choose the Non-Null Matrices 

Card 33.1 indicates that the load equation matrices are on magnetic file "LDTAPB** 
which begins with a DYLOFLEX header array (see sec. 6.4.1). Card 33.2 must be used if 
the matrices are on cards or if "LDTAPE” does not contain the DYLOFLEX header 
array. 

Card 33.1^Signal "LDTAPE” With DYLOFLEX Header Array 

Omit card 33.1 if the load equations are on cards or if "LDTAPE” does not contain 
the DYLOFLEX header array. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

DYLOFLX 

AlO 

Keyword indicates that "LDTAPE" begins with a DYLOFLEX header array. 

11-70 



Available for comments. 

■ L 


"LDTAP^” (card set 31.0) or "INTAPE” (card 4.4) must be defined in order to use 
card 33.1. When card 33.1 is used, all other BLOCK III cards (33.2 through 48.0) must 
be omifted. The DYLOFLEX header array and the load equations will be read from 
magnetic file "LDTAPE”. The DYLOFLEX header array at the beginning of "LDTAPE" 
specifies the sizes of the matrices and indicates which ones are not null and appear on 
the file (see sec. 6.4.1). 
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Card 33.2-Choose Non-Null Load Equation Matrices 
. Omit card 33.2 if card 33.1 was input. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

NMBl 

15 

- 

NMBl ,1*1,6 

If NMBl = 0 the corresponding load equations [M^] will not be 
read. It is assumed to be null. 

If NMBl = 1 the matrix [M^] will be read. 

6-10 

NMB2 

15 

11-15 

NMB3 

15 

16-20 

NMB4 

15 

21-25 

NMB5 

15 

26-30 

NMB6 

15 

31-35 

NC2B 

15 

If NC2B = 0 the load equation forcing function coefficient matrix 
ICgl or [C2] will not be read. It is assumed to 
be null . 

If NC28 = 1 the matrix or [C^] will be read. 

36-40 

NC3B 

15 

If NC3B = 0 the load equation convoluted forcing function coefficient 
matrix will not be read. It is assumed to be null. 

If NC3B = 1 the matrix or [C^] will be read. 

Note: If the number of panels (NPAN of Card 4.1) is > 0, then 

NC3B is set to zero. 

41-45 

NPHITB 

1 

1 

15 

If NPHITB = 0 the lifting panel contributions to the gust loads of 
designated gradual penetration load stations matrix 
[$] will not be read. It is assumed to be null. 

If NPHITB * 1 the matrix [$] will be read. 

Note: If the number of panels (NPAN of Card 4.1) * 0, then NPHITB 

is set to zero. 
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Card 33.3-Choose Non-Null Feedback Matrices 


Omit card 33.3 if card 33.1 was input or if IFDBAK = 0 on card 4.2. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

NSBl , i • 1 ,6 

If NSBi = 0 the load equation feedback coefficient matrix 

[S^] will not be read. It is assumed to be null. 

If NSBi = 1 the matrix [S^] will be read. 

1-5 

NSBl 

15 

6-10 

NSB2 

15 

11-15 

NSB3 

15 

16-20 

NSB4 

15 

21-25 

NSB5 

15 

26-30 

NSB6 

15 


Omit card sets 34.0 through 48.0 if the matrices are to be read from magnetic file 
"LDTAPE” (card set 31.0). 

Note: Format of matrices input on cards is detailed previously in section 6.3.1. 

Card Set 34.0 [M]]-Load Equation Matrix Corresponding to Generalized 

Coordinate Displacement 

Omit card set 34.0 if NMBl = 0 on card 33.2. [M] ] is of size NLD x NDOF. 

Card Set 35.0 [M 2 ] -Load Equation Matrix Corresponding to Generalized 

Coordinate Rate 

Omit card set 35.0 if NMB2 = 0 on card 33.2. [M2] is of size NLD x NDOF. 

Card Set 36.0 [M 3 ] -Load Equation Matrix Corresponding to Generalized 

Coordinate Acceleration 

Omit card set 36.0 if NMB3 = 0 on card 33.2. [M3] is of size NLD x NDOF. 

Repeat card set 37.0 through 48.0 NKVAL times; NK = 1, NKVAL. 

Card Set 37.0 [M 4 ]-Load Equation Matrix Corresponding to Convoluted 

Generalized Coordinate Displacement 

Omit card set 37.0 if NMB4 = 0 card 33.2. [M4] is of size NLD x NDOF. 

Card Set 38.0 [Ms]-Load Equation Matrix Corresponding to Convoluted 

Generalized Coordinate Rate 

Omit card set 38.0 if NMB5 = 0 on card 33.2. [M5] is of size NLD x NDOF. 
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Card Set 39.0 [M^J-Load Equation Matrix Corresponding to Convoluted 

Generalized Coordinate Acceleration 

Omit card set 39.0 if NMB6 = 0 on card 33.2. [Me] is of size NLD x NDOF. 

Omit card sets 40.0 through 45.0 if IFDBAK = 0 on card 4.2. 

Card Set 40.0 [Si]-Load Equation Feedback Coefficient Matrix for Generalized 
Coordinate Displacement 

Omit card set 40.0 if NSl = 0 on card 33.3. [Si] is of size NLD x NDOF. 

Card Set 41.0 [S 2 ]-Load Equation Feedback Coefficient Matrix for General 

Coordinate Rate 

Omit card set 41.0 if NS2 = 0 on card 33.3. [S2] is of size NLD x NDOF. 

Card Set 42.0 [S 3 ] -Load Equation Feedback Coefficient Matrix for General 

Coordinate Acceleration 

Omit card set 42.0 if NS3 = 0 on card 33.3. [S3] is of size NLD x NDOF. 

Card Set 43.0 [S 4 ]-Load Equation Feedback Coefficient Matrix for Convoluted 
Generalized Coordinate Displacement 

Omit card set 43.0 if NS4 = 0 on card 33.3. [S4] is of size NLD x NDOF. 

Card Set 44.0 [S 5 ]-Load Equation Feedback Coefficient Matrix for Convoluted 
Generalized Coordinate Rate 

Omit card set 44.0 if NS5 = 0 on card 33.3. [S5] is of size NLD x NDOF. 

Card Set 45.0 [Sej-Load Equation Feedback Coefficient Matrix for Convoluted 
Generalized Coordinate Acceleration 

Omit card set 45.0 if NS6 = 0 on card 33.3. [Sg] is of size NLD x NDOF. 

Card Set 46.0 {C 2 }-Load Equation Coefficient Matrix for the Forcing Function 

Omit card set 46.0 if NC2 = 0 on card 33.2. {C2} is of size NLD x 1. Input as a row. 

Card Set 47.0 {C 3 }-Load Equation Coefficient Matrix for the Convoluted Forcing 
Function 

Omit card set 47.0 if NC3 = 0 on card 33.2 or if NPAN > 0 on card 4.1. 

If there is a single forcing function (NFORC = 0 on card 4.1), {C3} is of size NLDx 1 
and is input as a row. 
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However, if the equations of motion are frequency dependent (NKVAL > 0 on card 4.1), 
{C3} is complex and must be input as a row of size 2*NLD with alternating real and 
imaginary parts. 

If there are multiple forcing functions (NFORC > 0 on card 4.1), [C3] is of size 
NLD X NFORC. 

However, if injaddition the equations of motion are frequency dependent (NKVAL > 0 
on card 4.1), [C3] is complex and must be input by rows with 2*NFORC numbers per 
row, alternating real and imaginary parts. 

Card Set 48.0 [$] -Gradual Penetration Load Equation Coefficient Matrix for the 
Convoluted Forcing Function 

Omit card set 48.0 if NPAN = 0 on card 4.1. [if>] is of size NLD x NPAN. 

[^] contains real numbers if NKVAL = 0 on card 4.1 and complex numbers if 
NKVAL > 0, indicating that the matrices are frequency dependent. When [<#>] is 
complex, it is input with alternating real and imaginary parts, 2* NPAN numbers to 
a row. 

Card Set 49.0-Termination of FINDRMS 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SEND 

AlO 

This card will cause the termination of the execution of FINDRMS. 




After this card is read, control of the L221 (TEV156) program will 




be returned to the main overlay, TEV156. 
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6^.4 PROGRAM SORTQLS 


All of the cards read by the second primary overlay SORTQLS are identified by a code 
word beginning in column one and have fixed field input of the required numbers. 

The recognized code words and possible parameters are shown in the following listing. 
The keyword parameters are underlined, and only the capital letter variables are 
relevant. 


Card set 

Keyword 

Variables 

51 

TITLE 

Miscellaneous titling information 

52 

TAPE 

Magnetic file '*NEWTPE” used for output 

53 

RESTArt 

Restart with MDOF, MFREQ, MSTATE, MPLTPE, 
MRTAPE, NGROUP, and NSETS 

54 

LOAD 

Load set number JSET 

55 

STATIC 

Static-elastic solution from load set JSET 

56 

PAIR 

Pair the standard and the static-elastic solutions 

57 

SEPARate 

Separate standard and the static-elastic solutions 

58 

SCALE 

Scale the axes with SCALX and SCALY 



Sort and retain the requested item: 



^ Q Generalized coordinates 



L Loads 

59 

SORT ^ 

S Spectra 



R RMS values 

60 

SORT 

A All items sorted and retained 

61 

$END 

Terminates SORTQLS execution 

Genera] information on 

the processing of the second primary overlay SORTQLS 


presented in the following paragraphs. 

Card sets 52.0 through 58.0 (TAPE, RESTA, LOAD, STATI, PAIR, SEPAR, and SCALE) 
simply set options or give directions to SORTQLS-no immediate action is taken. 
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The SORT card (card sets 59.0 and 60.0) causes SORTQLS to read the input magnetic 
file, sort the items specified, and write the magnetic file "NEWTPE”. Tlierefore, all 
desired options must be declared before the SORT card is input. 

SORTQLS will continue to process data cards and sort the specified items until it reads 
the $END card (card set 61.0) which will terminate the SORTQLS processing. 

The option cards LOAD, STATI, PAIR, SEPAR, and SCALE (card sets 54.0 
through 58.0) may be inserted any number of times. The last specified card set will 
remain in effect. 

The flow of SORTQLS card input is displayed in figure 6. 
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Card Set 50.0-SORTQLS Call Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

$S0RTqIs 

AlO 

This card initializes the call to the second primary overlay, 
SORTQLS. It introduces SORTQLS data cards. SORTQLS will perform 
the required sorting, and return control to the main overlay after 
Card Set 61 has been read. 

11-70 



Available for comments. 


Card Set 51.0-Title Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

TITLE 

AlO 

A keyword indicating a title card. 

11-70 

description 

6A10 

Any miscellaneous title description. 

The card will be read, printed on the output file, and ignored. 


Card Set 52.0-Tape Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

TAPE 

AlO 

A keyword indicating a tape definition card. 

11-20 

NEWTPE 

AlO 

The sorted plot vectors will be written in READTP/WRTETP format 
on the magnetic file "NEWTPE". 


Card set 52.0 should be input only once at the beginning of SORTQLS’s execution. 
Magnetic file "NEWTPE” will contain the sorted plot vectors. If this option is not used, 
then the default of "NEWTPE” is PLTPSD. 
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Card Set S3.0>Restart L221 (TEV156) 


Card set 53.0 is intended to restart L221 (TEV156) to sort a FINDRMS tape from a 
previous run. It is needed to specify information to SORTQLS which is normally passed 
through labeled common from FINDRMS. 

Card 53.1-Restart Card Variables 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

RESTArt 

AlO 

A keyword introducing a card of problem size variables. 

11-15 

MOOF 

15 

The number of degrees of freedom. 
(MOOF < 100) 

16-20 

MFREQ 

15 

The number of frequencies. 
(MFREQ >1) 

21-25 

MSTATE 

15 

The static-elastic solution. 

= 0 no static-elastic solution was previously found. 
= 1 a static-elastic solution was previously found. 

26-35 

MPLTPE 

AlO 

The magnetic file "IPLTPE" previously defined by FINDRMS. 

If MPLTPE is to be defined, then leave the variable blank, 

otherwise, the magnetic file name must be a left-justified, 1 to 7 
character file, name with the first character being alphabetic. 

36-45 

MRTAPE 

AlO 

The magnetic file which contains the load transfer functions. 

If MRTAPE is not to be defined, then leave the variable blank, 
otherwise, the magnetic file name must be a left-justified, 1 to 7 
character file name, with the first character being alphabetic. 

46-50 

NGROUP 

15 

The number of items contained on the magnetic file defined by MPLTPE 
for each frequency. 

(NGROUP > 3) 

51-55 

NSETS 

15 

The number of sets of load equations. 
(NSETS < 10) 

56-60 

ISET 

^15 

If NSETS = 1 the number of loads In load set 1. If NSETS Is 
greater than 1 this field Is to be left blank. 
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Read card 53.2 to specify the number of loads in each load set if NSETS of card 53.1 > 1. 
Card 53.2-Number of Sets of Load Equations 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

ISET^ 

1015 

ISET., 1=1, NSETS 

6-10 

ISETg 


where NSETS of Card 53.1 is the number of sets of load 




equations and must be £ 10. 

• 

• 

• 

• 


ISET^ contains the number of loads in the ith load set. 

46-50 



(ISET. < 100) 


Card Set 54.0-Load Set Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

LOAD 

AlO 

A keyword indicating that the Load Set to be processed follows. 

11-15 

JSET 

15 

The load set number. 




(JSET < NSETS) 


Card set 54.0 specifies the load set from which subsequent items will be sorted. If a 
LOAD card is not included, JSET = 1 is assumed. 

Card Set 55.0-Static Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

STATIC 

AlO 

A keyword requesting the processing of the static-elastic solution 
for a particular Load Set. 

11-15 

JSET 

15 

The load set number. 
(JSET < NSETS) 


Card set 55.0 specifies that subsequent items to be sorted will be from the static-elastic 
solution for the load set JSET. 
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Card Set 56.0-Pair Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

PAIR 

AlO 

A keyword requesting that standard and static-elastic responses 
be paired on the output file, "NEWTPE" 


: Card set 56.0 causes SORTQLS to sort the same specified items from both the standard 
‘ and static-elastic solutions (same load set) and write them on the magnetic file deHned 
by "NEWTPE” (card set 52.0), immediately following one another so they may easily be 
plotted on one frame. This feature is turned off by the SEPAR (card set 57.0) card. 

I 

i 

; Card Set 57.0-Separate Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SEPARate 

AlO 

A keyword requesting the separation of standard and static-elastic 




responses. 


Card set 57.0 will turn off card set 56.0 and will separate the standard and the 
static-elastic solutions. 

Card Set 58.0-Scale Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

r — - • 

SCALE 

AlO 

A keyword introducing SORTQLS Scale factors. 

ll-20 

SCALX 

ElO.O 

The frequencies will be multiplied by SCALX before being 
written on the magnetic file "NEVfTPE'*, (see Card Set 51) 
(SCALX > 0) 

21-30 

SCALY 

ElO.O 

The sorted Items will be multiplied by SCALY. 
(SCALY > 0) 
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Card Set 59.0-Sort Card 


Card 59. 1-Sort Q, L, S, or R 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SORT 

AlO 

A keyword requesting the sorting of item(s) from magnetic fne"INTAPE“ 

11-20 

a.k.s. 

or R 

AlO 

A variable indicating which item to sort from load set JSET. 
g^*s = sort the generalized coordinates, 

l^oads = sort the load transfer functions. 

Spectra = sort the output spectrum. 

^s = sort the root mean square values. 

21-25 

NSORT i 

15 

The number of items to be sorted. 

(NSORT <100) 

If NSORT = 0 all available items of the type specified above will be 
sorted (retained) and written on file "NEWTPE”. 


Card 59.1 is repeated for each variable which the user wishes to sort. 
Omit card 59.2 if NSORT of card 59.1 = 0. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

NUM^ 

1415 

NUM^ ; 1 = 1, NSORT 

6-10 

NUMg 


where NSORT is from Card 59.1. 


• 

• 


The specific numbers of items that are to be sorted (retained) as 

• 

• 


specified by the item type (Q, L, S, or R) on Card 59.1. 


n^^sort 


Repeat the card with 14 numbers per card until all elements are 




defined. 
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Card Set 60.0-^Sort 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SORT 

AlO 

A keyword requesting the sorting of items from file “INTAPE**. 

11-20 

MJL 

AlO 

This variable indicates what type of sorting will take place. 

Ait = will sort all A's, loads, and output spectra, from all load 




sets (both standard and static-elastic solutions if present). 


The SORTQLS magnetic output file will contain: 


Q’s for standard solution. 

Q’s for static-elastic solution if present. 

Loads for standard solution. 

Loads for static-elastic solution. 

Spectra for standard solution. 

Spectra for static-elastic solution. 

and, if IPLRMS = 1 (card 4.3 of FINDRMS) 


repeated 
NSETS times 

repeated 
NSETS times 


A values for standard solution. 

A values for static-elastic solution. 


repeated 
NSETS times 


Card Set 61.0-SORTQLS End Card 


qOLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SEND 

AlO 

A keyword which terminates the SORTQLS processing. 
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6.3.5 PROGRAM PLOTQLS 


All of the cards read by the third primary overlay PLOTQLS are identified by a code 
word beginning in colunin one and have fixed field input of the required numbers. 

The recognized code words and possible parameters are shown in the following listing in 
alphabetical order. The keyword parameters are underlined and only the capital letter 
variables are relevant. 


Card set 

Keyword 

Variables 

71 

DEFAUlt 

Use default grid limits 

70 

GRID 

Set grid limits to XMIN, XMAX, YMIN, and YMAX 

72 

LABEL 

Use miscellaneous plot labeling information 

75 

LINEAr 

Requests a linear-linear grid on the plot frame 

65 

LOAD 

Load set number JSET 

76 

LOG-Log 

Requests a log-log grid on the plot frame 

67 

PAIR 

Pairs the standard and static-elastic solutions 
Plot the requested item: 

f 

Q Plots the generalized coordinates 

77 

PLOT ^ 

L Plots the load transfer function 

S Plots the output spectrum 

^ R Plots the RMS load values 

78 

PLOT 

A Plots all items 

64 

RESTArt 

Restart with MDOF, MFREQ, MSTATE, MPLTPE, 
MRTAPE, NGROUP, and NSETS 

69 

SCALE 

Scales the axes with SCALX and SCALY 

68 

SEPARate 

Separates the standard and static-elastic solutions 

66 

STATIC 

Plot the static-elastic solution of load set JSET 

63 

TITLE 

Miscellaneous printout labeling information 

73 

X-AXIs 

Set x-axis label 
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74 Y-A Xls 


Set y-axis label 


79 SEND Terminate PLOTQLS processing 

General information on the processing of the third primary overlay PLOTQLS is 
presented in the following paragraphs. 

The cards DEFAUlt, GRID, LABEL, LOAD, PAIR, RESTArt, SCALE, SEPARate, 
STATIC, and TITLE simply set options or give directions to PLOTQLS-no immediate 
action is taken. 

The PLOT card (card sets 70.0 and 71.0) causes PLOTQLS to read the input magnetic 
file, sort the items specified, and prepare the plots. Therefore, all desired options must 
be declared before the PLOT card is input. 

PLOTQLS will continue to process data cards and sort the specified items until it reads 
the $END card (card set 79.0) which will terminate the PLOTQLS processing. 

The option cards DEFAUlt, GRID, LOAD, PAIR, SCALE, SEPARate, STATIc, and 
TITLE may be inserted any number of times. The last specified card set will remain 
in effect. 

The flow of PLOTQLS card input data is displayed in figure 7. 
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• Optional 


Figure 7.— (Continued} 














Card Set 62.0-PLOTQLS Unk Call Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

tPLOTqIs 

AlO 

This card initializes the call to the third primary overlay. 
PLOTQLS. It introduces PLOTQLS data cards. PLOTQLS will 
perform the required plotting and then will return control to 
the main overlay after Card Set 79 has been read. 


11-60 

description 

5A10 

Miscellaneous labeling information for plot identification. It 
should Include the user's name, phone number, nail stop, and organization. 


Card Set 63.0-Title Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

TITLE 

AlO 

A keyword indicating a TITLE card. 

11-70 

lescriptior 

6A10 

Any miscellaneous title description. 


Card set 63.0 will be read, printed on the output file, and ignored. 





Card Set 64.0-Restart L221 (TEV156) 

Card set 64.0 will restart L221 (TEV156) to sort a FINDRMS tape from a previous run. 
It is needed to specify information to PLOTQLS which is normally passed through 
labeled common from FINDRMS. 

Card 64.1-Restart Card Variables 


COLS. 

KEYWORD/ 

VARIABLE 

■ 

FORMAT 

DESCRIPTION 

1-10 

RESTArt 

AlO 

A keyword introducing a card of problem size variables. 

11-15 

MDOF 

15 

The number of degrees of freedom. 
(MDOF < 100) 

16-20 

MFREQ 

15 

The number of frequencies. 
(MFREQ >1) 

21-25 

MSTATE 

15 

The static-elastic solution. 

= 0 no static-elastic solution was previously found, 
= 1 a static-elastic solution was previously found. 

26-35 

MPLTPE 

AlO 

The magnetic file "IPLTPE" previously defined by FINDRMS. 

If MPLTPE is to be defined, then leave the variable blank, 

otherwise, the magnetic file name must be a left-justified, 1 to 7 
character file, name with the first character being alphabetic. 

36-45 

MRTAPE 

AlO 

The magnetic file which contains the load transfer functions. 

If MRTAPE is not to be defined, then leave the variable blank, 
otherwise, the magnetic file name must be a left- justified, 1 to 7 
character file name, with the first character being alphabetic. 

46-50 

NGROUP 

15 

The number of items contained on the magnetic file defined by 
MPLTPE for each frequency. 

(NGROUP > 3) 

51-55 

NSETS 

15 

The number of sets of load equations. 
(NSETS < 10) 

56-60 

I SET 

15 

If NSETS = 1, the number of loads in load set 1. If NSETS is 
greater than 1 then this field is to be left blank. 
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Read card 64.2 for the number of loads in each load set if NSETS of card 64.1 > 1. 


Card 64.2-Number of Loads in Each Load Set 


cots. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

ISET^ 

1015 

ISET., i « 1, NSETS 

6-10 

: 

• 

ISET2 

• 

■ 

• 


where NSETS of Card 72.1 is the number of sets of load equations 
and must be < 10. 

46-50 

ISETjO 


ISET^ contains the number of loads in the 1th load set. 
(ISET. < 100) 


Card Set 6S.0-Load Set Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

MO 

load 

AlO 

A keyword Indicating that the Load Set to be processed follows. 

11-15 

JSET 

15 

The load set number. 


' ■ 


(JSET S NSETS) 


Card set 65.0 specifies the load set from which subsequent items will be plotted. If a 
LOAD card is not included, JSET = 1 is assumed. 

Card Set 66.0-Static Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

MO 

STATIC 

AlO 

A keyword requesting the processing of the static-elastic solution 




for a particular load set. 

IM5 

JSET 

15 

The load set number. 




(JSET < NSETS) 


Card set 66.0 specifies that subsequent items to be plotted will be from the static*elastic 
solution of load set number JSET. 
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Card Set 67.0-Pair Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

PAIR 

AlO 

A keyword requesting that standard and static-elastic responses be 
paired on the plots. 


Card set 67.0 causes PLOTQLS to plot the same specified items from both the standard, 
and static-elastic solutions (of the same load set) on the same frame. This feature is 
negated by the SEPARate card (card set 68.0). 

Card Set 68.0-Separate Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SEPARate 

AlO 

A keyword requesting the separation of standard and static-elastic 
responses. 


Card set 68.0 will negate the PAIR option of card set 67.0, and will separate the 
standard and static-elastic solutions. 

Card Set 69.0-Scale Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SCALE 

AlO 

A keyword introducing PLOTQLS scalars. 

11-20 

SCALX 

ElO.O 

The frequencies will be multipl ied by SCALX before being plotted. . 
(SCALX > 0) 

21-30 

SCALY 

1 ElO.O 

1 

The plot items will be multiplied by SCALY before being plotted. 
(SCALY > 0) 
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Card Set 70.0-Grid Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

GRID 

AlO 

A keyword introducing PLOTQLS grid limits. 

U-20 

XMIN 

ElO.O 

The minimum x-axis grid limit. 

21-30 

XMAX 

ElO.O 

The maximum x-axis grid limit. 

31-40 

YMIN 

ElO.O 

The minimum y-axis grid limit. 

41-50 

YMAX 

ElO.O 

The maximum y-axis grid limit. 


Card set 70.0 establishes the limits of the plot card. If a GRID card is not input, the 
grid limits will be the lowest and highest frequencies and the minimum and maximum 
of the curves plotted on each frame. If a negative limit is specified, the program will 
- revert to the default for that particdlar limit. 

Card Set 71.0-Default Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

' l-io 

OEFAUIt 

AlO 

The keyword indicating that PLOTQLS will find its own grid limits 
for each frame. It also negates the GRID option of Card Set 70. 


Card Set 72.0-LabeI Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

l-io 

LABEL 

AlO 

A keyword introducing a plot label. 

' 11-50 

description 

4A10 

Any miscellaneous plot frame labeling information. 


The 40 label characters will be placed in the plot frame title block below the message(s) 
listed. Two label cards may be used. 


GEN COORD ^ 
LOAD TRANS > 
SPECTRUM 
RMS ^ 


NO. 


LOAD SET 


FROM THE STATIC-ELASTIC SOLUTION^ 
WITH THE STATIC-ELASTIC SOLUTION/ 


Note: The y-axis label is 

substituted for this 
line if input. 


40 characters of label card 1 
40 characters of label card 2 
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Card Set 73.0-X>Axis Label Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

X-AXIs 

AlO 

A keyword introducing an x-axis label. 

11-30 

description 

2A10 

Miscellaneous information to label the x-axis of the plot frame. 
The default x-axis label is FREQUENCY, 


Card Set 74.0- Y> Axis Label Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

Y-AXIs 

AlO 

A keyword introducing a y-axis label. 

11-50 

description 

4A10 

Miscellaneous information to label the y-axis of the plot frame. 


The y-axis label will serve as the first line in the title block. The default y-axis label 
and title line are: 


GEN COORD I 
LOAD TRANS > 
SPECTRUM 
RMS 


NO. 


LOAD SET 


Card Set 75.0-Linear Grid Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

LINEAr 

AlO 

A keyword requesting a linear - linear grid on the plot frame. 
Note: This is the default grid type. 


Card Set 76.0-Log Grid Card 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

LOG-Log 

AlO 

A keyword requesting a log - log grid on the plot frame. 
Note: The default grid type is linear - linear. 
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Card Set 77.0-Plot Card 


Card 77.1-Plot Q, L, S, or R 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

PLOT 

AlO 

A keyword requesting the plotting of item(s) from magnetic f11e''INTAP! 

11-20 

a.L.s. 

or R 

AlO 

A keyword indicating which item to plot from load set JSET. 
g[*s = plot the generalized coordinates. 

U)ads = plot the load transfer functions. 

Spectra = plot the output spectra. 

^s » plot the root mean square values. 

21-25 

NPLOT 

15 

The number of items to be plotted. 

(NPLOT S 100) 

If NPLOT = 0 all available items of the type specified above 
win be plotted. 


Card 77.1 is repeated for each variable the user wishes to plot. 
Omit card 77.2 if NPLOT of card 77.1 = 0. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

NUM^ 

1415 

NUM. ; 1 = 1. NPLOT 

6-10 

NUM^ 


where NPLOT is from Card 77.1. 




The specific numbers of items that are to be plotted as specified by 

* 

• 


the item type (Q, L, S, or R) on Card 77.1. 


^UHnpLOT ^ 


Repeat the card with 14 numbers per card until all elements are 




defined. 



Card Set 78.0-Plot A 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

PLOT 

AlO 

A keyword requesting the plotting of items from magnetic file "INTAPE* . 

11-20 

Ate 

1 

AlO 

This keyword requests that all Q's, loads, and output spectra, 
from all load sets (both standard and static-elastic solutions if 
present). 


The plots will be in the following order: 

Q’s for standard solution. 

Q’s for static-elastic solution if present. 

Loads for standard solution. 

Loads for static-elastic solution. 

Spectra for standard solution. 

Spectra for static-elastic solution. 

and, if IPLRMS = 1 (card 4.3 of FINDRMS) 

A values for standard solution. 

A values for static-elastic solution. 

Card Set 79.0-PLOTQLS End Card 


repeated 
NSETS times 

repeated 
NSETS times 


repeated 
NSETS times 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 


1-10 



$end 

1 

AlO 

A keyword which terminates the PLOTQLS processing. 


Card Set 80.0-Termination of L221 (TEV156) 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 


1-10 

$QUIT 

AlO 

A keyword which terminates the L221 (TEV156) program. 
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Requrrements 
or Function 

Keywords and/or Variables 

j Card 
! Format 

Reference 
Card Set (CS) 


i 

SFREQuency 






I 

AlO 

1.1 


sflTLe 






A 10 

1.2 


SCORE 

HAXCOR 





; Aio 

1.3 


SCHECk 






A 10 

1.4 


SFIWDnas 






AlO 

2.0 




BLOCK i - 

- Options and Constant Data 


i 

1 



LABEL 

7A10 

! 3.0 

Options 

NOOF 

NFREQ 

NPAN 

NKVAL NFORC 

NDOFD 


615 

4.1 

Options 

IFREQ 

ISPEC 

ISTATE 

IFDBAK I DAMP 

INTZRO 


615 

4.2 

Output Options 

IPRINT 

IRMSPR 

I PUNCH 

ICKPRT I PL RMS 

LPOOF 


615 

4.3 

I/O Files 

INTAPE 

IPLTPE 

IRTAPE 

IFTAPE 



4A10 

4.4 

Wagner Function 

A, 

Bi 

Cl 

Di Ql 

Bi 

Yi 

7E10 

5.1 

Kiissner Function 

Ai 

B 2 

C 2 

D 2 0L2 

62 

Y 2 

7E10 

5.2 

Input Const. 

z 

VEL 

L 

C TABFS 

TABS 

FOSCAL 

7E10 

6.0 

ISTATE = 2; CS 4.2 

Columns to 

be zeroed out for static- 

•elastic solution 



1415 

7.1 

NDOFD > 0; CS 4.1 

Columns and rows of Equations of Motion and columns of 
Load Equations to be zeroed out 

1415 

7.2 

ISPEC » 2; CS 4.2 

K 

n 





2E10 

8.1 

ISPEC < 0; CS 4.2 
IFTAPE undefined 
CS 4.4 

Tabular input spectrum frequencies 

7E10 

8.2 


6J5.S SUMMARY OF CARD INPUT DATA 
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Requirements 
or Function 

Keywords and/or Variables 

Card 

Format 

Reference 
Card Set (CS) 

ISPEC < 0; CS 4.2 
NFORC = 1; CS 4.1 
IFTAPE undefined; 
CS 4.4 

Single tabular input spectrum 

7E10 

8.3 

ISPEC < 0; CS 4.2 
NFORC ^ 2; CS 4.1 
IFTAPE undefined 
CS 4.4 

Multiple tabular input spectrum 

7E10 

8.4 

IDAMP = 0 
CS 4.2 

9sD 

lElO 

9.1 

IDAMP > 0 
CS 4.2 

Structural damping for each degree of freedom 

7E10 

9.2 

IFDBAK = 1 
CS 4.2 

Standard Feedback Coefficients 

7E10 

10.0 - 10.5 

IFDBAK = 2 
CS 4.2 

Special Feedback Coefficients 

7E10 

^ 10.6 - 10.8 

IFREQ « -1 
CS 4.2 

Frequencies (cps) 

7E10 

n.i 

IFREQ = 0 
CS 4.2 

Frequencies (rad/sec) 

7E10 

11.2 

IFREQ > 0 

CS 4.2 1 

First frequency and an increment (cps) 

2E10 

11.3 

LPDOF > 0 i 

CS 4.3 

Frequencies where loads/degree of freedom is printed 

7E10 

1 11.4 



X- 


Requirements 
or Function 

1 

Keywords and/or Variables 

Card 

Format 

Reference 
Card Set (CS) 


SEQUAtion 

AlO 

12.0 

If tape from 
DYLOFLX system 

DYLQFIx 

AlO 

13.1 

If matrices on cards, 
or INTAPE not from 
DYLOFLX 

Choose non-null equations of motion matrices 

1015 

13.2 

IFDBAK > 1 
CS 4.2 

Choose non-null feedback matrices 

615 

13.3 

Matrices on cards 

Equations of motion matrices and associated frequencies (if any) or 
equations of motion feedback matrices 

7E10 

14.0 - 30.0 


S 


II-In*tructions and Matrices to Calculate the General Coordinates Frequency Respoaahs 
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Requi rements 
or Function 

Key Words and/or Variables 

Card 

Format 

Reference 
Card Set (CS) 


$L0AD NLD ICORR T LDTAPE 

A10,2I5, 

E10,A10 

31.0 

ICORR > 0; CS 31.0 
NFORC = 1; CS 4.1 

Load number. Other load numbers to correlate with first item 

1115 

32.0 

If tape from 
DYLOFLX system 

DYLOF.Ix 

AlO i 

33.1 

If matrices on cards, 
or LDTAPE not from 
DYLOFLX 

Choose non-null load equation matrices 

1015 

33.2 

IFDBAK > 1 
CS 4.2 

Choose non-null feedback matrices 

615 

33.3 

Matrices on cards 

Load equation matrices or load equation feedback matrices 

AlO 

34.0 - 38.0 


|END 




BLOCK Ill-Instructions and Matrices to Calculate the Loads Frequency Responses 



Requirements 
or Function 


Optional 


Keywords and/or Variables 


Reference 
Card Set (CS) 


SSORTqIs 


TITLE - Title description 


Tape card 


Restart capability if generalized coordinates and load frequency response tape 
previously saved 


load JSET 


STATI C JSET 


£A1R 


SEPAR ate 


SCALE - Scale card A5,5X; ElO.O; E10.0 


SORT - Items and number of items to sort 


Specific numbers of items to be sorted 


50.0 


52. 


55.0 


57.0 

58.0 




ST 


I 

1 

1 


NSORT > 0 
CS 59.1 





















s 


Requi rements 
or Functions 

Keywords and/or Variables 

Reference 
Card Set (CS) 


SPLOTqls User’s nane/phone muiber/mail stop/orgainzation 

62.0 


TITLE - Title description 

63.0 


Restart capability if generalized coordinates and load frequency response tape 
previously saved 

64.0 


Specific items and controls for plotting 

65.0 - 78.0 


SEND 

79.0 





$QUIT - Terminates the L221 (TEV156) program 

80.0 


PLOTQLS-Writes a Special Magnetic File for Plotting on COMpSO 



6.4 MAGNETIC FILES INPUT DATA 


The input matriceB to the L221 (TEV156) program can be obtained from three magnetic 
(tape or disk) flies. These flies, "INTAPE”, "LDTAPE", and "IFTAPE", are discussed in 
sections 6.4.1, 6.4.2, and 6.4.3, respectively. 

The use of these magnetic flies to contain the input matrices will depend on whether or 
not they have been deflned through the use of either the tape option card (see card 4.4 
of sec. 6.3) or the load instruction card (see card set 31.0 of sec. 6.3). If any of the 
magnetic flies have not been deflned, then the associated input matrices will be read 
from input cards. 

All input magnetic files are in the READTP/WRTETP format.* 


6.4.1 ’TNTAPE” 

By default FINDRMS will read the equations of motion matrices in card sets 14.0 
through 30.0 (sec. 6.3). If, however, the user specifies a file name in the variable 
INTAPE on card 4.4, the matrices will read from magnetic file "INTAPE” (see fig. 8) in 
the READTP/WRTETP format. 

FINDRMS rewinds "INTAPE” before attempting to read any matrices. Several input 
parameters determine which matrices are contained on "INTAPE". These parameters 
are NDOF, NPAN, and KNVAL of card 4.1; IFDBAK of card 4.1 and the non-null 
matrix indicators in cards 13.2 and 13.3. 

For the matrices appearing in figure 8: 

{HEADER} Will be present only if the keyword DYLOFLEX was read in 

card 13.1 

[S] ],[S2],...,[Ss] Will not be present if IFDBAK = 0 on card 4.2 

{fji} and [$] Will replace {Cg} if NPAN > 0 on card 4.1 

{FREQM} and the frequency-dependent matrices ([M4], [M5], [Mg], {C3}, or [5*]) will 
not be present if NKVAL = 0. 

The equations of motion are followed on "INTAPE” by an end-of-flle. If LDTAPE in card 
set 31.0 is deflned with the same name as INTAPE, FINDRMS will skip the end-of-file 
and read the load equations from the second file of "INTAPE”. Additional load sets may 
appear in logical files three, four, etc. 


Clemmons, R. E. : Programming Specifications for Modules of the Dynamic Loads Analysis 
System to Interface With FLEXSTAB. NASA contract NASl-13918; BCS-G0701. 
September 1975. (Internal Document.) 
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Card 13.1 contains DYIOFLX 


{HEADER} 


If NKVAL > 1 
{FREQH} 


(HEADER) is 


[M.] 

is 

{NDOF X NDOF) 

[s]] 

is 

{NDOF X NDOF) 

{FREQM) 

is 

(NKVALx 1) 

{C^} 

is 

(NDOF X 1) 


or 

(NDOF X NFORC) 

{C3}* 

is 

(NDOF X 1) 


or 

(NDOF X NFORC) 

fi>} 

is 

(NPAN X 1) 

[$? 

i s 

(NDOF X NPAN) 

* Complex 

if NKVAL > 1 


If IFDBAK>0 
[S,].[S2],[S3], 
[S.].[Ss].[S,] 


If NPAN > 0 

{f£}. [^] 


Notes: 

• NDOF.NPAN, NKVAL, NFORC 
are obtained from 
Card 4.1. 

• IFDBAK is obtained 
from Card 4.2. 

• If NFORC >1 , {Cilis 
null and {f^} and 
[J] are not read. 


If NKVAL > 1 


If NPAN <0 

1 . ' ,*f 


^ I i=2, NKVAL 


£nd-of-File 


Figure 8.— Contents of the File "INTAPE" 







The contents of the DYLOFLEX header matrix are described in the following listing. 


Position Format 


Description 

1 

AlO 

Contains the DYLOFLEX header keyword, 7HDYLOFLX 

2 

AlO 

Program name/version number of the program that created the 
tape 

3 

AlO 

Date of the run that created the tape, lOH yr/mo/da 

4 

no 

Contains the number of degrees of freedom, NDOF 

5 

no 

Contains the number of load equations, NLD 

6 

no 

Contains the number of gradual penetration panels, NPAN 

7 

no 

Contains the number of k values used for the multiple 
frequency-dependent matrices option, NKVAL 

8 


Not used 


9 

EIO.O 

Contains the velocity (true air speed), VEL 

10-20 


Not used 


21 

no 

NULMAT 

(1); if > 0, [Mj] will be read 

22 

no 

NULMAT 

(2); if > 0, [M2] will be read 

23 

no 

NULMAT 

(3); if > 0, [M3] will be read 

24 

no 

NULMAT 

(4); if > 0, [M4] will be read 

25 

no 

NULMAT 

(5); if > 0, [M5] will be read 

26 

no 

NULMAT 

(6); if > 0, [Mg] will be read 

27 

no 

NULMAT 

(7); if > 0, [C2] will be read 

28 

no 

NULMAT 

(8); if > 0, [C3] will be read 

29 

no 

NULMAT 

(9); if > 0, [f£] will be read 

30 

no 

NULMAT (10); if > 0, [5] will be read 


Note: The program will check the values given for NDOF, NLD, NPAN, NKVAL, and 
VEL in the DYLOFLEX header matrix against the same variables read from 
cards. All values not compatible will cause a fatal error, except VEL which 
causes a v/arning error. 
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6.4.2 ’’LDTAPE” 


By default, FINDRMS will read the load equations in card sets 34.0 through 48.0 
(sec. 6.3). If, however, the user specifies a file name in the variable LDTAPE on card 
set 31.0, the matrices corresponding to those card sets will be read from magnetic file 
"LDTAPE^’ (see fig. 9) in READTP/WRTETP format. 

FINDRMS will rewind '^LDTAPE’* if LDTAPE does not have the same name as INTAPE 
(card 4.4) or LDTAPE of the previous load set. Otherwise, FINDRMS will skip one 
end-of-file mark and try to read the set of load equation matrices from the next logical 
file on "LDTAPE”. 

Several input parameters determine which matrices are contained on LDTAPE. These 
parameters are NDOF, NPAN, and NKVAL of card 4.1; IFDBAK of card 4.2; and the 
non-null matrix indicators in cards 33.2 and 33.3, 

For the matrices appearing in figure 9: 

(header [ Will be present only if the DYLOFLEX keyword was encountered 

in card 33.1, The contents are similar to those appearing on 
INTAPE described in section 6.4.1 

[Si ],[S2] . ...[Sg] Will not be present if IFDBAK = 0 

[5] Will replace {C3} if NPAN > 0 

The frequency-dependent matrices ([M4], [M5], [Mg], {C3}, or [^]) will not be present 
if NKVAL = 0. 

6.4.3 *TFTAPE” 

When ISPEC < 0 on card 4.2, FINDRMS will read the tabular input spectrum(s) on 
card 8.3 or 8.4 (sec. 6.3). If, however, the user also specifies a file name in the variable 
IFTAPE on card 4.4, the tabular input spectrum(s) will be read from "IFTAPE” 
(see fig. lO) in the READTP/WRTETP format. 

FINDRMS will rewind "IFTAPE”and then read one matrix for each spectrum required. 
A spectrum is needed for each combination of (I,J) where 1=1, NFORC and J = I, 
NFORC. Because of symmetry, only the upper right triangular combinations are input. 
Input row 1 with J = 1, NFORC; then row 2 with J = 2, NFORC; then row 3 with J = 3, 
NFORC ; ...; and finally row NFORC with J = NFORC. 

Each spectrum contains NSP elements, but because the elements are complex, the array 
must be of size 2NSP x 1. 


no 




Card 33.1 
contains DYLOFLX 


{HEADER} 




[ 1 ] 

1 II 1 

tn rs> 

1 II 1 





{HEADER} is 

(30 

X 1) 

[W,] 

is 

(NLD 

X NDOF) 


is 

(NLD 

X NDOF) 

{(Tl) 

is 

(NLD 

X 1) 

^ -k 

or 

(NLD 

X NFORC) 

{C,} 

is 

(NLD 

X 1) 


or 

(NLD 

X NFORC) 

w* 

is 

(NLD 

X NPAN) 


If IFDBAK>0 



♦Complex if NKVAL>0 

Notes: 

• NOOF.NPAN.NKVAL.NFORC 
are obtained from 
Card 4.1 . 

• IFDBAK is obtained 
from Card 4.2. 

• If NFOROl , rtz) 
is null and [?] 
is not read. 


If NPAN <0 If NPAN >0 



I Repeat 
^i=2,NKVAL 


Repeat for each load set, 
including End-of-File. 


Figure 9.— Contents of the File "LOT APE" 
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Size 


Matrix 


NSP X 1 
2NSP X 1 
2NSP X 1 

2NSP X 1 

2NSP X 1 
2NSP X 1 



where 

NSP = IlSPECi = the number of 
frequencies defining » 
the input spectrum 

NFORC = the number of 

excitation forces 




.. .-4 


Figure 10.— Contents of the File "IFTAPE" 



6.5 OUTPUT DATA 


6k6.1 PRINTED OUTPUT DATA 

The input data for each of the L221 (TEV156) modules will be echo-printed except 

where noted; i.e., printing of FINDRMS input matrices is optional. Other printed output 

matrices will be discussed by module. 

Program FINDRMS 

liie RMS values and zero crossings will be printed for each FINDRMS solution. 

Optionally, the following can also be printed. 

• Generalized coordinates, load transfer functions, and output spectrum at each 
frequency (see option IPRINT of card 4.3, sec. 6.3) 

• Correlation between different loads (see option ICORR of card set 30.0) 

• Input matrices; load equations and equations of motion (see option ICKPRT of 
card 4.3) 

• Matrices of intermediate calculations for program checkout (see option ICKPRT of 
card 4.3) 

• RMS load values. A, for the last IRMSPR frequencies (see option IRMSPR of 
card 4.3) 


Program SORTQLS 

After each set of card input data, SORTQLS will print one line for each item it writes 
on magnetic file "NEWTPE”. The messages are intended to show the program’s progress 
if a failure occurs. 


Program PLOTQLS 

After each set of card input data, PLOTQLS will print a one line message for each item 
plotted. The messages are intended to show the program’s progress if a failure occurs. 

•.5.2 MAGNETIC FILES OUTPUT DATA 

FINDRMS will write on as many as three magnetic files; "IPTLPE”, "IRTAPE” (figs. 11 
and 12, respectively), and "SCRATCH”. SORTQLS, however, will only write on one 
Magnetic file, "NEWTPE” (fig. 13). PLOTQLS writes on magnetic file "TAPE99”. 

All matrices contained on these magnetic files are written by the WRTETP subroutine.^ 


'Clemmons, R. E.; Programming Specifications for Modules of the Dynamic Loads Analysis 
System to Interface With FLEXSTAB. NASA Contract NASl-13918; BCS-G0701, 
September 1975. (Internal Document.) 
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"IPLTPE” 


The magnetic file "IPLTPE” is written when it is requested on card 4.4 of FINDRMS. 
"IPLTPE”, if defined, will contain the arrays calculated in FINDRMS (fig. 11). 


Size 

Matrix 

(NFREQ 

X 1) 

{FREQ} 

(NDOF 

X 1) 

{QMAGR} 

(NLD 

X 1) 

{QMAGL} 

(NLD 

X 1) 

{SPEC}* 

(NDOF 

X 1) 

{QMAGR} 

(NLD 

X 1) 

{QMAGL} 

(NLD 

X 1) 

{SPEC}* 

• 

(NDOF 

X 1) 

• 

• 

{QMAGR} 

(NLD 

X 1) 

{QMAGL} 

(NLD 

X 1) 

{SPEC}* 


End-of-File 



{FREQ} is in cycles per second 


for 

FREQ(l) 


for 

FREQ(2) 


for 

FREQ(NFREQ) 

Each Load set generates two 
solutions if the static elastic 
solution was requested with one 
file created for each solution 


If IPLRMS = 1 on Card 4.3, {A} which contains the A values 
will be written on "IPLTPE" immediately after {SPEC} for each 

frequency. 


Figure 1 1 .—Contents of the File "IPLTPE” 



I 


"IRTAPE” 

The magnetic file "IRTAPE” is written when it is requested on card 4.4 of FINDRMS. 
"IRTAPE”, if defined, will contain the arrays calculated in FINDRMS (fig. 12). 


(NFREQ 

X 1) 

{FREQ} 


(NDOF 

X 1) 

{Qr} 

a 

(NDOF 

X 1) 

{qT} 


(NLO 

X 1) 

{SOHPp} 


(NLO 

X 1) 

{SOWj} 

- 

(NDOF 

X 1) 

{Qr} 

a 

(NDOF 

X 1) 

{Qt} 


(NLD 

X 1) 

{^r} 


(NLD 

X 1) 

{SDflCj} 

• 

m 

(NDOF 

X 1) 

• 

{Qr} 

- 

(NDOF 

X 1) 

{Qt> 


(NLD 

X 1) 

{ISOHCr} 


(NLD 

X 1) 

{ §UMC j } 

a 



End-of-F11e 







Each Load set generates two 
solutions If the static elastic 
solution was requested with one 
file created for each solution 


for 

FREQ(l) 


for 

FREQ(2) 


for 

FREQ(NFREQ) 


{FREQ} Is in cycles per second 

{Qp} and {Qj) are the real and 
Imaginary parts of {Q} 


(SUMCJ and (SUMCj) are the real 
and Imaginary parts of {5DRC}. 


Figure 12.— Contents of the File "IRTAPE" 


115 


I 


III I III! 


I II 


II I I I I I 


I I III! II II I 


1 1 HIM II I iiii 


"NEWTPE” 

The magnetic file "NEWTPE” is written whenever SORTQLS is run. "NEWTPE” will, 
contain alternating matrices of independ€ht and dependent variables to be plotted by 
another program (see fig. 13). 


Size 

Matrix 

(NFREQ X 1) 

{FREQ} 

(NFREQ X 1) 

{ITEMi} 

(NFREQ X 1) 

{FREQ} 

(NFREQ X 1) 

{ITEM 2 } 

• 

• 

(NFREQ X 1) 

• 

{FREQ} 

(NFREQ X 1) 

{ITEM^} 


End-of-File 


{FREQ} = SCALX * {FREQ} 

{ITEM^.} = SCALY times the ith item sorted 
by SORTQLS. 


Figure 13.-Contents of the File "NEWTPE" 



6.5.3 PUNCHED CARDS 


There is an option in FINDRMS that will cause the program to punch on cards the RMS 
load values, {A}, and number of zero crossings, {Nq}, for each solution. This option is 
requested by IPUNCH of card 4.3 (sec. 6.3). 

The matrices {A} and {Nq} will be punched on cards with the format of 7E10.3 with 
each matrix starting at the beginning of a card (flg. 14). 



Figure 14.— Format of jAj and Punched on Cards 


H7 



6,6 RESTRICTIONS 


General: 

• Maximum core size allowed for L221 (TEV156) is 200000g. 

• All input variables must be within the ranges defined in the card input data 
(sec. 6.3) 

FINDRMS: 

• A maximum of 10 loads may be correlated with any particular load 

• Loads per degree of freedom may be printed at a maximum of 10 different 
frequencies 

• A maximum of 10 load sets in one data case 

Use of the multiple forcing function option (app. C) results in the following restrictions: 

• No gradual penetration 

• No feedback 

s No loads correlation 

• No static-elastic solution 

• Generalized coordinates and load frequency response functions may be plotted for 
only the first excitation force 

• {^ 2 } and {C 2 } are null 

• Loads per degree of freedom may not be printed 

• Kiissner lift growth function must = 1 

• Only the tabular input spectrum can be used 
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6.7 DIAGNOSTICS 


All errors detected by L221 (TEV156) will result in the printing of a diagnostic error 
message. These messages are self-explanatory and are of the following formats: 

******** FATAL ERROR n 

DIAGNOSED WHILE EXECUTING ROUTINE name 
other lines describing the error 

******** WARNING ERROR n 

DIAGNOSED WHILE EXECUTING ROUTINE name 
other lines describing the error 

where n is the error number (from 1 to 19), and name is the name of the routine in 
execution when the error was detected. 

All errors are fatal unless they have been diagnosed as warnings. The following list is a 
brief description of each error. 

1. Premature end-of-flle was encountered in the card input file. 

2. Keyword is not appropriate. 

3. Keyword is not recognizable to L221 (TEV156). 

4. Tested parameter is not within its specified range or it is in conflict with other options 
(see the card input data, sec. 6.3). 

5. Illegal magnetic file definition. 

6. Required core is found to be too large. 

7. READTP error (see sec. 6.7.1). 

8. WRTETP error (see sec. 6.7.2). 

9. PRINTM/PRNTCM error (see sec. 6.7.3). 

10. Improper matrix size. 

11. FSF error. The returned value i contains the number of file marks remaining to be 
skipped when an end-of-information was encountered. If n files were to be spaced, 
then (n-i) files were spaced before the error was detected. 

12. [SUMM] matrix was found to be singular when routine CGLESM was called; CGLESM 
solves AX = B for X. 

13. Extrapolation is required over the tabular input spectra. 

14. VARDIM routine errors (see sec. 6.7.4). 
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15. FETEDIT errors (see sec. 6.7.5). 

16. Rudder limiting (RLIMIT) did not converge. This is nonfatal error and execution 
will continue (see app. B). 

17. Unable to sort or plot the RMS load values, A, because they were not contained in 
file 'TPLTPE”. This is a nonfatal error and execution will continue. 

18. Number of requested plots is incorrect. This is a nonfatal error and execution will 
continue. 

19. Unable to plot 0. on log grid. The current plot will be omitted. This is a nonfatal 
error and execution will continue. 


6.7.1 READTP ERROR CODES 


Error code 


0 No errors are detected during reading. 

lOOO+I A forward space file (FSF) error occurred, where I is the 
number of file marks remaining to be skipped when an 
end-of-information was encountered. 

2 Number of matrices or files to be skipped, before reading 
starts, is less than zero. 

3 Dimensioned number of rows in the matrix is less than 
zero. 

3000+1 A forward space record (FSR) error occurred, where I is 
the number of records remaining to be skipped when either 
an end-of-file or end-of information was encountered. 

4 Number of rows in the matrix is greater than the 
dimensioned row size in the program. 

5 Name check failed. 

6 Number of rows in the matrix (M) times the number of 
columns (N) is greater than the buffer size; or M * N ^ 0. 

7 End-of-file was read. If it occurs while reading the matrix 
ID, no information is stored in the user's area. If it occurs 
while reading the matrix, the ID information will be stored. 
Note that the records will always be in pairs and an 
end-of-file should always be encountered with the ID record. 
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6.7.2 WRTETP ERROR CODES 


Error code 


0 No errors are detected during writing. 

lOOO+I FSF error occurred, where I is the number of file marks 
remaining to be skipped when an end-of-information was 
encountered. 

2 Number of matrices or Ales to be skipped, before writing 
starts, is less than zero. 

3 Dimensioned number of rows in the matrix is less than or 
equal to zero. 

3000+1 FSR error occurred, where I is the number of records 
remaining to be skipped when either an end of file or 
end-of-information was encountered. 

4 Actual number of rows in the matrix is greater than the 
dimensioned number of rows in the matrix. 

6 Number of rows in the matrix (M) times the number of 
columns (N) is greater than the buffer size. 


6.7.3 PRINTM/PRNTCM ERROR CODES 


Error code 


0 No error was detected. 

3 Number of dimensioned rows is less than or equal to zero. 

4 Actual number of rows is greater than the number of 
dimensioned rows. 

6 Actual number of rows times the actual number of 
columns is less than or equal to zero. 


6.7.4 VARDIM ROUTINE ERROR CODES 

DELETE error code = 0 Successful. 

= -2 Name of the array is not contained within the 

VARDIM library. 

INITIR error code ^ 0 Successful. 

= -1 Previous array of the same name was destroyed before 

the new one was established. 
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= 1 Maximum number of catalog entries was exceeded. 

= 2 One of the array dimensions is zero. 

= 3 Blank common storage was exceeded. 

= 4 Routine STARTR was not called beforehand. 

LOCATE error code = 0 Successful. 

= -1 STARTR was not called beforehand. 

= -2 Name of the array is not contained within the 

VARDIM library. 

6.7.5 FETEDIT ERROR CODES 
FETADD error code = 0 Successful. 

= -1 File has been previously defined, no action taken. 

= 1 An illegal tape number. 

= 2 Length of either the buffer or array was too small. 

= 3 Maximum number of files has already been defined 

(49 allowed). 

FETDEL error code = 0 Successful. 

= -1 File was not defined, no action taken. 

= 1 An illegal tape number. 
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7.0 SAMPLE PROBLEM 


The sample problem shown in this section is a vertical gust analysis of a t3rpical large 
subsonic jet transport airplane. The theoretical math model consists of two rigid body 
degrees of freedom (vertical translation and pitch) and six wing elastic normal modes. 
Unsteady aerodynamics was modeled only for the wing with the Doublet Lattice 
method. The rigid body aerodynamics representation (rigid body stability derivatives) 
was obtained from corrected wind tunnel data. 

The load equations were generated using the force summation method. 
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Sample Problem Card Input Data 


smcQufticv 






1.0 

0ATA4 

2 

trtuc 

00U8LET 

LATTICE UNSTEADY AEAO CHECK CASE 


1.1 

0A7A4 

3 

tmtE 







l.l 

DATA4 

4 

STIUf 


HULTIFLE FKEQUENCV 

OEFENOENT 

MATRICES 


1.1 

0ATA4 

7 

$T|U6 







1.1 

DATA 4 

8 

STITtE 

THE EQUATIONS OF 

MOTION WILL 

BE BEAD FROM THE MAGNETIC FILE 

1.1 

0ATA4 

9 

SMTLE 

EOMTFE. 






l.l 

0ATA4 

10 

»T|TL6 







l.L 

0ATA4 

11 

•TITLE 

THE LOADS EQUATIONS WILL BE 

READ FROM 

THE MAGNETIC FILE 

l.l 

0ATA4 

12 

•TITLE 

LOOTFE. 






l.l 

DA TA4 

13 

•TITLE 







1.1 

DATA4 

14 

•EINOMS 







2.0 

DATA4 

15 

FIIEQUENCV 

OEFENOENT AEAOe 8 

DEG. OF FKEEDOMt 8 LOADS# 210 

FREQUENCIES 

3.0 

DATA4 

16 

8 710 0 

20 0 

0 




4.1 

0ATA4 

17 

0 

t 0 

0 0 

1 




4.2 

0ATA4 

18 

1 

0 0 

1 0 

0 




4.3 

EMPOS 

1 

EONTFE 

FLOTTFE 






4.4 

EHFOS 

2 

ta 

0. 

0. 

0, 

0. 

0. 

0. 

WAGNER 

0ATA4 

21 

1. 

0. 

0. 

0. 

0.* 

0. 

0. 

KUSSNER 

0ATA4 

22 

laO 

829.5332 

2500. 





6.0 

0ATA4 

23 

OaO 







9.1 

DATA4 

24 

• 05 

• 1 

• 2 

• 3 

•4 

• 45 

•5 

W-1 11 

20ATA4 

25 

•55 

.4 

• 65 

• 7 

• 8 

• 9 

1. 

H-2 11 

20ATA4 

26 

l.l 

1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

H-3 ll 

20ATA4 

27 

l«8 

1.9 

2.0 

2.1 

2.2 

2.4 

2.7 

W-4 11 

20ATA4 

28 

3*0 

3.3 

3.6 

3.9 

4.2 

4.5 

4.8 

W-5 11 

2DATA4 

29 

5.0 

5.25 

5.5 

5.75 

6.0 

6.25 

6.50 

W-6 11 

20ATA4 

30 

6.75 

7.0 

7.25 

7.50 

7.75 

6.0 

8.2$ 

W-7 11 

20ATA4 

31 

8* 50 

8.75 

9.0 

9.25 

9.50 

9.75 

10. 0 

W*8 11 

20ATA4 

32 

10.25 

10.50 

10.75 

ll.O 

11.25 

11.50 

11.75 

W-9 ll 

20ATA4 

33 

12.0 

12.25 

12.50 

12.75 

13.0 

13.29 

13.50 

H-lOtl 

20ATA4 

34 

13.75 

14.0 

14.25 

14.50 

14.75 

15.0 

15.25 

W-llll 

20ATA4 

35 

15.50 

15.75 

16.0 

16.25 

16.50 

16.75 

17.0 

H>1211 

20ATA4 

36 

17.25 

17.50 

17.75 

18.0 

18.25 

18.50 

18.75 

W-1311 

20ATA4 

37 

19.0 

19.25 

19.50 

19.75 

20.0 

20.25 

20.5 

W-14U 

20ATA4 

38 

20.75 

21.0 

21.25 

21.5 

21.75 

22,0 

22,25 

W-1511 

2DATA4 

39 

22.5 

22.75 

23.0 

23.25 

23.50 

23.75 

24.0 

H*161L 

20ATA4 

40 

24.25 

24.50 

24.75 

25.0 

25.25 

25.50 

25.75 

H-1711 

20ATA4 

41 

26.0 

26.25 

26.50 

26.75 

27.0 

27.25 

27.50 

W->18U 

20ATA4 

42 

27.75 

28.0 

28.25 

28.50 

28.75 

29.0 

29.25 

H-I9U 

20ATA4 

43 

29.50 

29.75 

30.00 

30.25 

30.50 

30.75 

31.0 

W-20 1 1 

20ATA4 

44 

31.25 

31.50 

31.75 

32.0 

32.25 

32.50 

32.75 

W-2UI 

20ATA4 

45 

33.0 

33.25 

33.50 

33.75 

34.0 

34.25 

34.50 

W-2211 

2DATA4 

46 

34.75 

35. 

35.5 

36. 

36.5 

37. 

37.5 

W-23U 

20ATA4 

47 

38. 

38.5 

39. 

39.5 

40. 

40.5 

41. 

W-2411 

20ATA4 

48 

41.5 

42. 

42.5 

43. 

43.5 

44. 

44.5 

W-25U 

20ATA4 

49 

45. 

45.5 

46. 

46.5 

47. 

47.5 

48. 

W-26U 

20ATA4 

50 

48.5 

49. 

49.5 

SO. 

50.5 

51. 

51.5 

M*2711 

20ATA4 

51 

52. 

52.5 

93. 

53.5 

54. 

54.S 

55. 

W-2811 

20ATA4 

52 

55.5 

56. 

96.5 

57. 

57.5 

58. 

58.9 

W*2911 

2DATA4 

53 

42.5 

69 . 

T5. 

•2. 

88 . 

94. 

99. 

M-3011 

20ATA6 

54 

OEOMOrfOH 






12.0 

DATA4 

55 

1 

0 1 

1 1 

0 0 

1 



U.l 

0ATA4 

56 

•L080 

• 

0 

tOOTFC 




31.0 

0ATA4 

57 

0 

0 1 

1 1 

0 0 

1 



33.1 

0ATA4 

58 

•END 







49.0 

DATA4 

59 
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Card Input Data (Concluded) 


SSOKTOLS 



50*0 

EMPOS 

3 

TITLE 



51*0 

EHPOS 

4 

TITLE 

SAMPLE 

SORTQLS EXECUTION* 

S1*0 

EMPOS 

5 

TITLE 



Si.O 

EMPOS 

6 

SORT 

0 

1 

59.1 

EMPOS 

7 

3 



59*2 

EMPOS 

8 

SORT 

LOADS 

t 

59.1 

EMPOS 

9 

1 



59.2 

EMPOS 

10 

iCNO SORTQLS EXECUTION 


6U0 

EMPOS 

11 

IPLOTOLS 

USERS NAME / 

telephone no* / MAIL STOP / ORG* 

4i2*0 

EMPOS 

12 

title 



63*0 

EMPOS 

13 

TITLE 

SAMPLE 

PLOTQLS EXECUTION. 

63.0 

EMPOS 

14 

title 



63*0 

EMPOS 

IS 

PLOT 

0 

1 

77.1 

EMPOS 

16 

3 



77.2 

EMPOS 

17 

LABEL 

EXAMPLE PLOriNG LABEL 

72*0 

EMPOS 

16 

PLOT 

LOADS 

1 

77.1 

EMPOS 

19 

1 



77.2 

EMPOS 

20 

SEND PLOTQLS EXECUTION 


79.0 

EMPOS 

21 

SQUIT THE 

TEV156 PPOGRAMMING SYSTEM 

ao.o 

oataa 

60 
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Sample Problem Printed Output 


• PROGRAH t22l VERSION El 

• 8 EG INNING EXECUTION ON THE 8CS 6600a 

• DATE OF RUN IS 76/09/29* 

• TIME OF RUN IS 2UU.l3a 


UFREQUENCY 1.0 
ISTIUE DOUBLET LATTICE UNSTEADY AERO CHECK CASE !•! 
ISTITLE X.l 
ISTITLE multiple frequency dependent MATRICES i.l 
liUTLE 1.1 
(STITLE the equations OF MOTION WILL BE READ FROM THE MAGNETIC FILE. 1.1 
ISTITLE EOMTPE. 1.1 
IsriTLC lal 
(STITLE the LOADS EQUATIONS HILL BE READ FROM THE MAGNETIC FILE 1.1 
isiitle LODTPE. 1.1 
ISTITLE 1.1 
(SFINORMS 2.0 
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ffNOtMS* THC PflOGRAN TO FIND TH€ PSD RMS AND NUMftCR 
OF ZFRO CROSS INGSf IS NOW IN EXECUTION* 

ACCUMULATED CP SECONDS • .9A50 


THE INITIAL CARD INPUT TO FINDRHS FOLLOWS. 

DATA CASE LABEL ICARO SET 3. OK 

I FREQUENCY DEPENDENT AERO» B DEG* OF FREEOOMt S LOADSt 210 FREQUENCIES I 

PROBLEM SI2E ICARO A*l)« 

NUMBER OF DEGREES OF FREEDOM INOOFI • 8 


NUMBER 

OF 

FREQUENCIES 

INFRSQI « 

210 

NUMBER 

OF 

K- VALUES 

INXVALI • 

20 


FREQUENCY DEPENDENT MATRICES WILL BE READ* 

INTERPOLATION WILL BE USED TO FIND 
THE MATRICES AT THE SOLUTION FREQUENCY* 

OPTIONAL CAPABILITIES ICARO 4*21* 

FREQUENCIES (IFREO) WILL BE READ FROM CARDS IN RAO./SEC* 

SPECTRUM NUMBER IISPECI I WILL BE USED* 

THE SAME DAMPING FACTOR WILL BE USED FOR EACH FREEDOM (lOAHP *GT*OI* 

integration WILL BEGIN AT FREQIHt AND NOT ZERO I INT2R0*EQ* 1). 

OUTPUT OPTIONS (CARO 4*31* 

THE GEN* COORD. AND LO* TRANSFER FOR EACH FREQUENCY AS WELL AS THE 
FINAL RMS AND NO. OF ZERO CROSSINGS WILL BE PRINTED C IPRINT *EQ* 1) • 

PRINT ONLY THE FINAL RHS VALUES 1 1 RMSPR.EQ.O J • 

THE INPUT MATRICES WILL BE PRINTED I ICKPRT •GT.OI* 

INPUT TAPE OPTIONS ICARD 4*4). 

THE INPUT MATRICES WILL BE READ FROM EORTPE 
THE PLOT DATA KILL RE WRITTEN ON PLOTTPE 

WAGNER FUNCTION COEFFICIENTS ICARD 5* 1) • 

.lOOOEtOl 0* 0. 0* 

KUSSNER FUNCTION COEFFICIENTS ICARD 5.2). 

*1000E«>01 0* 0* 0* 


0 * 0 * 0 * 

0 * 0 * 0 * 


THE CORE REQUIRED FOR TA&ULAR INPUT SPECTRUPIS) IS 00001177^6 


INPUT CONSTANTS (CARD SET 6l« 
Z - *100QE4>01 

VEL • *829SE^03 

i. « •2500E^C4 


The STRUCTURAL 

DAMPING FACTOR (CARO 9« 11 

USED FOR ALL DEGREES OF 

FREEDOM - 

0. 

The frequency hatrix (freq) 

IN CYCLES/SECONO. 




.7956E-02 

♦1592E-01 

«3l83E-0l 

•4775E-01 

•63666*01 

•7162E*01 

•79586-01 

•0754E-01 

.9549E-01 

• 103584^00 

.11146400 

.12736400 

•14326400 

•1592E400 

.1751E^00 

. 1910E^OO 

.2069E>00 

•2228E400 

.2387E400 

•2546E400 

.27066400 

*2865E^00 

.3024E>00 

.318364-00 

•33426400 

•35016400 

•33206400 

•42976400 

• 4775E4.00 

.5252Ef00 

.5730E4.00 

.62076400 

•66856400 

•71626400 

.76396400 

.7958E+00 

.8356E^00 

.8754E400 

.91516400 

.95496400 

.99476400 

.10356401 

• 1074»E^01 

• lIl4Et01 

• 11S4E4^01 

.11946401 

.12336401 

•12736401 

.13136401 

. 1353E^01 

. 1393E^01 

•1432E401 

.14726401 

•15126401 

.15526^01 

.15926401 

•1631E«01 

. 167lE^0l 

.I7UE401 

.17516401 

. 17906401 

.18306401 

.18706401 

• 1910E4-01 

• l950Et0l 

.1989E401 

.20296401 

•20696401 

.21096401 

•21496401 

• 21P8E«^01 

«2228Ef01 

• 226884^01 

.2308E401 

.23486401 

•23878401 

•24276401 

*2467E^0l 

• 2507E^0l 

• 2546E4-01 

.2S66E40L 

•26266401 

•26666401 

•2706E401 

.2745E4-01 

.2785£^Ol 

• 282SE4-01 

.286SE401 

.29056401 

.29446401 

.2984E401 

•302^E^01 

• 3064E^01 

.3l04£^0l 

•31436401 

.31836401 

.32236401 

.32636401 

•3302E*01 

«33^2E4-01 

• 3382E4'01 

•34226401 

•34626401 

•35016401 

.35416401 

.3501E1-OI 

• 3621E + 01 

•3661E^0L 

•3700E401 

•37406401 

.37806401 

•382064Q1 

• 3a60Et>01 

• 3099E4.OI 

• 3939t4-0l 

•3979E401 

.40196401 

.40586401 

•40986401 

•^13BE^0l 

. 41788^01 

«4216E401 

.42576401 

.42976401 

.43376401 

.43776401 

•A417E^01 

.A456E*0l 

.449664.01 

.45366401 

•45766401 

•46156401 

.46556401 

•469SE>01 

.47356^0l 

.4775E401 

•48146401 

.48546401 

.4894E401 

•49346401 

.4974F*0I 

• 50136^01 

.50536401 

.50936401 

•5133E401 

.51736401 

•52126401 

.52528^01 

.5292F^0l 

.5332E401 

.53716401 

.54116401 

•5451E401 

•54916401 

•S53IEf0l 

.5570E*01 

.56506401 

.57306401 

.58096401 

.58896401 

•59666401 

.6048E4-01 

«6127E*01 

.62076401 

.62876401 

.63666 fOl 

•64466401 

.^5256401 

•6605E401 

• 6685E4'01 

.67646401 

•68446401 

•69236401 

.70036401 

.70826401 

.7l62E*Ol 

•7242E^0l 

.73216401 

.74016401 

.74806401 

.75606401 

.76396401 

•7719E*01 

• 7799E*01 

.78786401 

.79506401 

•80376401 

.81176401 

•B196E401 

«8276E«01 

•83868^01 

•84356401 

•8515E401 

.85946401 

.86746401 

•87546401 

•8833E^0l 

• 89l3£t0l 

.69926401 

•90726401 

.91516401 

.92316401 

•93LIE401 

•9947Ef01 

• 1098E^02 

.1194E402 

•13056402 

•14016402 

•1496 6402 

•1576E402 
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l•f•\MTION 12.0 I 

t 10111001 13.1 I 

CHOOSE MATRICES TO t£ READ FOR OE-N. COORD. iCARO SET 131 »T CODE OF 0 OR 1. 
MATRICES ASSUMED NULL AND NOT READ IF CODE >0. 



m 

H3 Hh H9 

Mb 

SI 

S2 

$3 $4 

S5 $4 

1 

0 

1 1 1 

0 

0 

0 

0 0 

0 0 

C2 

C3 

FtL PHIT 






0 

1 

0 0 






THE 

CORE 

REQUIREO TO 

PlHO 

THE 

GEN. 

COOPO. 

IS 0000121S43 


MATtlX — — — — 

LABa Ml 

Sf2E 8 ROWS 

AND 

8 COLUMNS 






COLUMN 

1 

2 

3 

4 

5 

4 

7 

8 

ROW 1 

0. 

0. 

0. 

0. 

Or 

Or 

Or 

0. 

ROW 2 

0. 

0« 

0. 

Or 

Or 

Or 

Or 

Or 

ROW 3 

0. 

0. 

0* 

0. 

Or 


Or 

Or 

ROW 4 

0. 

0* 

0. 

•9483E«03 

Or 

Or 

Or 

Or 

ROW 5 

0. 

0. 

Or 

0. 

•4589E«>04 

Or 

Or 

Or 

ROW 6 

0. 

0. 

Or 

Or 

Or 

•3413E404 

Or 

Or 

ROW 7 

0. 

0. 

Or 

Or 

Or 

0. 

r6143E^0^ 

Or 

ROW 8 

0. 

0. 

Or 

Or 

Or 

Or 

Or 

• 9772C^04 
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HATHIX 

LABEL H3 

SIZE 8 ROWS AND 8 COLUMNS 


COLUMN 


1 

2 

3 


4 

3 

6 

7 

8 

ROW 

1 

•9788E^03 

0. 

0* 


0* 

0. 

0* 

0. 

0« 

ROW 

Z 

0« 

• 1729E4-09 

0* 


0. 

0. 

0. 

0 . 

0. 

ROW 

3 

0. 

0* 

•2396F#02 

0. 

0* 

0* 

0. 

0. 

ROW 

A 

0. 

0. 

0, 


.2517E+02 

0. 

0. 

0« 

0* 

ROW 

5 

0, 

0. 

0. 


0. 

•9367E4'01 0. 

0« 

0* 

ROW 

6 

0* 

0. 

0. 


0. 

0« 

•8987E^02 

0. 

0. 

ROW 

7 

0. 

0. 

0. 


0. 

0* 

0* 

• 6493E4'01 

0. 

ROW 

MATRIX 

LABEL 

8 

0* 

0« 

0. 


0* 

0. 

0. 

0. 

.89756^0l 


FREOH 









SIZE 


1 ROWS AND 

20 COLUMNS 








COLUMN 

1 

2 

3 

4 

$ 

6 

7 

6 

9 

10 


«4886E-01 

•aoooE+00 

•1698E^01 

•3594E^01 

•6241E«01 

•7988E^01 

.1098E>02 

•1348E^02 

• 1599E4-02 

•L897£«02 


.2197E+02 

•2596E+02 

.3I95E^02 

•3894E+02 

•4294Et02 

•4593E+02 

•5l92Ef02 

• 6690E4'02 

• 7489E4>02 

•998S£^02 


THE FREQUENCY DEPENDENT FREQUENCY MATRIX IFREOMI IN CYCLES/SECOND. 


•7776E-02 

.1273E*'00 

.2702E>00 

.5720E>00 

.9933E^00 

• 1271E^01 

•1748E>01 

•2146E>01 

.2544E+01 

.30l9£^0l 

•3496E*0l 

• 4I32E«>01 

•5085E*0l 

•6198E>01 

•6833E>01 

• 7310E'^0l 

•6264E>0i 

. 10976^02 

.1192E^02 

.IS69E>02 



Note: The equations of motion input matrices (M4, M5, and C3) have been omitted 
from this document. 


FREQUENCY DEPENDENT AERO. 8 DEG. OF FREEDCM* 8 LOADS* 210 FREQUENCIES 


Note: The frequency-dependent equations of motion input matrices (M4, M5, and C3) 
have been omitted from this document. 


FREQC 11 - .7S98E-02 CYCLES/SEC. ( .5000E-01 RAO./SEC.I 


RESPONSE 


FAGNITOCE 


PHASE ANGLE 


•10T13E*0? 

-•91767E-03 

.2636AE-02 

•15313E-03 

•88676E-(H 

.l35blE-0«» 

.89028E-0A 

-.27395E-OA 


.58339E*02 

29277E-0A 

-•1688AC-01 

-•707AAE-03 

-•AA332E-03 

-•2858AE-03 

-.3A107E-O3 

.108S6E-03 


.5931SE»02 

.91813E-03 

.17089E-01 

.72A02E-03 

.6A9A0E-03 

.2Bb\bE-03 

•35327E-03 

•1L196E-03 


• 13B92E«>01 

• 3173SE«'01 
.A8673E»0I 
«99235C^01 
•48994 E^Ol 
.47S98E»0l 
.49671£>0I 
•18L80E>01 


FREO( 2) - .lS92E-0i CYCLES/SEC* t .lOQOE^OO RAO./SEC.I 


RESPONSE 


MAGNITUDE PHASE ANGLE 


.17351E«01 

-.92075E-03 

.39171E-02 

•21431E-03 

•134S7E-03 

.33663E-04 

.U663E-03 

-•35889E-04 


•29034E902 

•32373E-04 

-.341S8E-01 

-.14381E-02 

-•13023E-02 

-•S7364E-03 

-.69594E-03 

.22082E-03 


.29086E.402 

.92L32E-03 

•34422E-01 

.14540E-02 

.13094E-02 

.57663E-03 

.70565E-03 

.22372E*03 


•15111E901 

•31064E401 

•48264E401 

•48603E401 

.48193E401 

•47708E401 

•48784£«0l 

•17319E«0l 



Note: Frequencies 3 through 208, including the associated frequency-dependent 

equations of motion input matrices (M 4 , M 5 , and C 3 ), have been omitted from 
this document. 


FREQI2C9) * .M96E>C2 CYCLES/SEC. ( •9400E>02 RAO./SEC*) 


RESPONSE 


HAGNtnOE PHASE ANGLE 


1 

2 

3 

4 
3 
6 

7 

8 


.74946E-03 
•49438E-06 
•3445fle-02 
•24296E-02 
• 57369E-03 
•220O1E-O3 
.63770E-02 
-•32289E-02 


•27962E-03 

•16401E-06 

•11457E-02 

•22180E-02 

♦19767E-04 

18003E-03 

.39910E-02 

.938356-03 


.799936-03 

.520886-06 

.363136-02 

.328976-02 

.574036-03 

.29052E-C3 

.752296-02 

.33625e-02 


.357106^00 
.320326^00 
.32100E4-00 
.73990£^00 
.34442 E-Ol 
.56146E>01 
.55920£^00 
•28588E^0l 


PR6Q(210) - .1S76E>02 CYCLES/SEC. ( .9900Et02 RAO./SEC.) 


RESPONSE 


MAGNITUDE PHASE ANCLE 


1 

2 

3 

4 

5 

6 

7 

8 


.681126-03 

•4S448E-QG 

.314736-02 

.228056-02 

.568776-03 

•19282E-03 

.57701t-02 

-.27317E-02 


•26C99E-03 

•14746E-Q6 

.11022E-02 

.203006-02 

.479416-04 

12941E-03 

.35661E-02 

.93460E-03 


.729416-03 
.477816-06 
.33347E-02 
•305316-02 
.570786-03 
.232226-03 
.67832E-02 
•28871 E-02 


.365926^00 

.31375£^00 

•336856^00 

• 72736E«>00 
•840908-01 
•569216^01 
.55357EfOO 

• 26 tl9E«'01 



mOAO 8 0 LOOTPE 3I.0 I 

LOTAPE WHICH CCNTAINS THE LOAD EQUATiQNS IS DEFINED AS LOOTPE 

The file defined by intape ieohtpe i will be replaced by lotape. 

THE NUMBER OF LUAC EQUATIONS IN SET I IS 0 


I 0 0 1 1 I 0 0 1 33.1 1 

CHOOSE MATRICES TO BE READ FOR LOADS SQLUTIOH tCARD SET 331, 

HATRICES ASSUMED NULL AND NOT READ IF CODE % 0. 


HBl 

H62 

MB3 NBA 

PBS HB6 

SBl 

SB2 

SB3 

SBA 

SB3 

SB6 

0 

0 

1 1 

1 0 

0 

C 

0 

0 

0 

0 

C28 

C3B 

PHITB 









0 10 


CORE REQUIRED TO FIND LOADS AND RMS VALUES IS 0000121601 


MATRIX 


LABCL 

SIZE 

HB3 

8 ROWS AND 

8 COLUMNS 






COLUMN 


1 

2 

3 

4 

5 

6 

7 

8 

ROW 

1 

•2020E«02 

.81986+04 

.12476+02 

.79256+01 

.46946+01 

.26026+01 

-.65856+00 

-.15456+01 

ROW 

2 

•1638E^03 

.3571E+05 

•46426+02 

-.5952E+01 

-. 12836+02 

-.22796+01 

•1503E+01 

.27266+01 

ROW 

3 

.A709E+03 

.92576+04 

.3624E+02 

-.95306+01 

16636+02 

.54726+02 

.9569F+01 

-.10936+02 

ROW 

A 

.25906-02 

.14566+01 

.25906-02 

.25906-02 

.25906-02 

.25906-02 

•2590E-02 

.25906-02 

ROW 

5 

•2590E-C2 

.88176+00 

•11826-02 

.48816-03 

-.67006-04 

-.2612E-03 

-.54006-03 

-.19166-03 

ROW 

6 

.25906-02 

.17416+00 

.69796-04 

-.30176-04 

-.29176-03 

.65026-03 

.31996-03 

-.28066-03 

ROW 

7 

.25906-02 

-.4028E+00 

-.1499E-03 

.58306-04 

.7573E-04 

-.28526-03 

-.43206-04 

.49496-04 

ROW 

8 

.25906-02 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


Note: The load equations input matrices (M4, M5, and C3) have been omitted from this 

document. 


mcQumcr oc^ehoeut 6 oeg. cf frieccf, a FucQuEiictES 


RESUtTS ffm lOAO SET ttO« t 


Note: The frequency-dependent load equations input matrices (M4, M5, and C3) have 
been omitted from this document. 


FREQI 1) • •79S0E-C2 CVCUES/SEC. ( .SCOOE-Oi RAO./SEC.I 

VALUE OF INPUT SPECTRUM • •81993E^C3 0. (LEN./SEC.I*P2/CPA0«PLEN«I 


LOADS MAGNITUDE PHASE ANGLE OUTPUT SPECTRUM 

IRAQ*) L0A0/(RA0./LEN.) 


1 -.32099E*01 *IA992E^02 

2 -•L0633E^02 .312$3E>02 

3 -•l8671Et02 •37249E^02 

4 -.^60516-04 -.37756E-03 

5 -•67357E-04 -.37767E-03 

6 -•689 73E-04 -•37 7 77E-03 

7 -•70295E-04 -.37782E-03 

8 -•69372E-04 -.37778E-03 


*15332E^02 *178l7E^0i 
•3304IE»02 •189ese^01 
•41667E>02 *20354E^0l 
•38329E-03 .45392E^0l 
«38363E-03 •4S3S9E^01 
• 36401E-03 «453l8E4>0i 
•38430E-03 .45284E^0l 
•36410E-03 «45308E»0l 


•19274Et06 

•69514E^06 

*1423SE^07 

•12046E-03 

•12067E-03 

.I2091E-03 

*l2109E-03 

.12097E-03 


FREOC 21 « .13928-01 CVCLES/SEC. i .lOOOE^OO RAO./SEC.I 

value OF INPUT SPECTRUM « «86555E^03 0« C LEN./SEC. |4*2/l RAO.*LEN« } 

LOADS MAGNITUDE PHASE ANGLE OUTPUT SPECTRUM 

CRAO.I LOAO/I R AO./L EN« I 


I 


-.43l79E^0l 

••12936E^02 

-•21940E«C2 

-.3U53E-04 

-•38871E-04 

••43343E-04 

->• 488448-04 

-•44943E-04 


• 30491E«>02 
•63927E^02 
•76721E^02 
-.75153E-03 
-.75192E-03 
-.732076-03 
-•75197E-03 
-.75204E-03 


•30756E^02 
•65226E«02 
.797966^02 
.73220E-03 
.752836-03 
.753326-03 
•73334E-03 
.75339E -03 


.17117E401 

.17708EOI 

.18493E^0l 

•46703E^Ot 

*46634E^0l 

•46548E»01 

«46476E^01 

•46527E«>01 


•81873E404 
•3b623E«07 
•35113E^07 
.489736-03 
.490S5E-03 
•49120E-03 
•49140 E-03 
.49128E-03 


8 



Note: Frequencies 3 through 208, ^clu^ng the ^sociated frequency-dependent load 
equations input matrices (M 4 , M 5 , and C 3 ), have been omitted from this 
document. 


FfiE0(209) « •U96E^C2 CYCLES/SCC. I «9900E>02 RAD«/SEC.) 


VALUE OF INPUT SPECTRUM • «10676E«>00 0. 


(LEN*/$£C.)P«2/IRA0«PLEN«i 


LOADS 


MAGNITUDE 


PHASE ANGLE 
(RAO . i 


OUTPUT SPECTRUM 
LOAO/CRAD./LEN.) 


1 •89A68E^02 

2 -.ISTlbE^Oa 

3 -•A9232E+02 
A -.2A838E>C0 

5 -.A1540E-01 

6 -.A5255E-01 

7 -.80A31E-02 

8 -.17153E-01 


-.32225E»02 

-.ieaiSE^03 

•2985AE>03 

-.19A65E^OO 

-.89769E-02 

-.148376-01 

-.419446-02 

-.83998E-02 


.950956*02 

.24515E*03 

.3C258E*03 

.315576*00 

.424996-01 

.47563E-01 

.936376-02 

. 183086-01 


.593756*01 

.401656*01 

.173426*01 

.380636*01 

.335446*01 

.345446*01 

.367916*01 

.349876*01 


.965416*03 

•64161E*04 

.977386*04 

.106316-01 

.192826-03 

.241516-03 

.936036-05 

.357846-04 


FR60(210) « .1576E*02 CVCLES/SEC. ( .99006*02 RAO./SEC.) 


VALUE OF INPUT SPECTRUM « .979236-01 0. 


(LEN./SEC. ) 4*2/1 RAO.*LEN.» 


LOADS 


MAGNITUDE 


PHASE ANGLE 
(RAO.I 


OUTPUT SPECTRUM 
LOAD/(RAD./LEN.) 


1 .900986*02 

2 -.169486*03 

3 -.125856*02 

4 -.258046*00 

5 -.421986-01 

6 -.447536-01 

7 -.829046-02 

8 -.172926-01 


-.313156*02 

-.178306*03 

•308416*03 

-.200446*00 

-.100406-01 

-.146616-01 

-.492456-02 

-.662576-02 


•95405E*02 

.245996*03 

.308676*03 

.326746*00 

.433766-01 

.471006-01 

.964276-02 

.185186-01 


.594616*01 

.395236*01 

•161166*01 

.380206*01 

.337526*01 

.345666*01 

.367766*01 

.350756*01 


•691306*03 

.592566*04 

.932976*04 

.10454E-01 

.18424E-03 

.217236-03 

.910506-05 

.335786-04 



FREQUENCY DEPENDENT AERO, 


8 DEG. CF FREECON* 8 LOADS* 210 FREQUENCIES 


RESULTS FOR LOAD SET NO. 1 


RMS ZERO CROSSINGS 

(PER LENGTH) 


.53591E+03 
.l5bl2E»04 
• 16769E»0A 
•12932E*00 
.A4672E-01 
•IT^SeE-Ol 
.1229<JE-01 
.12010E-01 


.10644E-02 

.U629E-02 

•22369E-02 

.6S833E-02 

•3381SE-02 

.63952E-02 

.189^0E-02 

•2S276E-02 


tiENO 


49.0 1 


FlNORNS IS FINISNEO. 

ACCUNULATEO CF SECONDS • 23.7090 



I ISOftTOtS 


50*0 


I 


$OKTQtS» THC PMOtthn TO SORT PLOT VECTORS FftON FINORMS IS NOW IN EXECUTION. 
ACCUH. CF SEC. • 23.7090 


$1.0 ) 

51.0 ) 

51.0 ) 

59.1 ) 


CORE REQUIRED TO SORT IS 0000045602 

SORT GEN. COORD 3 FROM STANDARD SOLUTION OF LOAD SET 1 


ititle 

ITlTtE SAMPLE SQRTQtS EXECUTION. 

ITITLE 

(SORT Q 1 

THE NUMBER OF ITEMS TO BE SORTED (NSQRT) • I 
THE NUMBERS THAT ARE TO BE SORTED. 

3 


(SORT LOADS 1 59.1 ) 

THE NUMBER OF ITEMS TO BE SORTED (NSQRT) • 1 

THE NUMBERS THAT ARE TO BE SORTED. 

1 

CORE REQUIRED TO SORT IS 0000045602 

SORT LOAD TRANS t FROM STANDARD SOLUTION OF LOAD SET 1 


(SEND SORTQLS EXECUTION 6U0 ) 

SORTQLS IS FINISHED. 

ACCUM. CP SEC. • 25. 3560 



CiPLOTQLS USERS |«A^4E / TELEPHONE NO. / NATL STOP / OR.G. 


62.0 


I 


PLOTQLS. the program USED TO PLOT DATA VECTORS 
IS NON IN EXECUTION. 

ACCUMULATED CP SECCNOS • 29.3600 


63.0 ) 

63.0 ) 

63.0 ) 

77-1 ) 


THE CORE REQUIRED TO PLOT 1 S 00000703AQ 


PLOT GEN. COORD 3 FROM STANDARD SOLUTtON OF LOAD SET I 


I TITLE 

ITITLE sample PLOTQLS EXECUTION. 

ititle 

I PLOT 0 I 

THE number OF ITEMS (NPLOT) TO BE PLOTTED - 1 

THE NUMBERS THAT ARE TO BE PLOTTED. 

3 


72.0 J 

77.1 ) 

the core required TO PLOT IS 0000070340 

PLOT LOAD TRANS 1 FROM STANDARD SOLUTION GF LOAD SET 1 


ILABEL example PLOTING LABEL 

IPLOT LOADS 1 

THE NUMBER OF ITEMS INPLOT) TO BE PLOTTED - I 
THE NUMBERS THAT ARE TO BE PLOTTED. 

1 


I SEND PLOTQLS EXECUTION 79.0 I 

PLOTQLS IS FINISHED. 

ACCUMULATED CP SECONDS - 27.9190 




tiQUlT THE TEV156 fROGRAHNlNG SYSTEM 


80*0 


I 


8L0TQI.S tS ftHtSHED* 

ACCUMULATED C^ SECONDS • 27.9S00 
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LOAD SET 



FREQUENCY t CPS 1 

OEM. COOffO. NO. 3. LOAD SET 1 




APPENDIX A 


FREQUENCY ARRAY FORMULATION 


FINDRMS can automatically prepare five different distributions of frequencies 

FREQi; i = 1. NFREQ 

over which the power spectral density (PSD) solution will be found. The option ISPEC 
(card 4.2) chooses one of the five distributions shown below. 

Given: Initial frequency, FREQ( 1 ) or cj | (card 11.3) 

Frequency increment, DELTAF or Aco (card 1 1 .3) 

Number of frequencies, NFREQ or n (card 4.1) 


If ISPEC=1 

U)j 

= CO, 

+ jAco 

^n-4 

= ^n-5 

+ 200Aco 

^^n-3 

= ^n-4 

+ 200Aco 

^n-2 

= ‘^n-3 

+ 1 OOOAco 

"n-1 

= ^^n-2 

+ 2000ACO 

^n 

= ^n-1 

+ 3000Aco 

If ISPEC=2 

COj 

= CO, 

+ jAco 

OJj 

= COj., 

+ 2Aco 

^^n-4 

= "n-5 

+ 200Aco 

‘^n-3 

= "n-4 

+ 200Aco 


= "n-3 

+ 1 OOOAco 

^n-\ 

= "n-2 

+ 2000ACO 

cOn 

= "n-1 

+ 3000ACO 


where j=l. ..., (n-5) 


2 

where j=l, ..., k, and k= - (n-5), truncated 
where j=k+l, ..., (n-5) 
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IfISPEC=3 


COj 

= wj 

+ jAco 

where j=l, ..., k, and k=^(n-5), truncated 

COj 

= COj.i 

+ 2Aco 

3 

where j=k+l, .... k 2 , and k 2 =-(n- 5 ), truncated 

COj 

= «j.i 

+ 4Aco 

where j=k 2 +l, ..., (n-5) 

^^11-4 

= ‘^n-S 

+ 200Au) 


‘^n-S 

= ‘‘^n-4 

+ 200ACO 


^n -2 

= ^n-3 

+ lOOOAco 


^n -1 

= ‘^n-2 

+ 2000ACO 



= ‘^n-I 

+ 3000AOJ 


IfISPEC=4 

^j 

= OJi 

+ jAco 

where j=l, k, and k=^(n-15), truncated 


= 

+ 2Aco 

3 

where j=k+l, k 2 , and k2=~(n-15), truncated 

COj 

= Wj -1 

+ 4Aco 

where j=k 2 +l, n-15 

COj 

= COj.i 

+ 6 Aco 

where j=n- 14, n-5 

^^n-4 

= '^n-5 

+ 2 OOA 0 J 


‘^n-S 

= ^n-4 

+ 2 OOA 0 } 


“n -2 

= ^n-3 

+ lOOOAw 


‘^n-l 

= “n -2 

+ 2000Aco 


"n 

= "n -1 

+ 3 OOOA 0 J 
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If IFSPEC=5 



=: 

CO, 


,5AclJ 


= 

UJj., 

+ 

Acj 

"j 


Wj.l 

+ 

lAco 

«j 


Wj., 


3 

<I 


= 

COj.| 

+ 

6Aco 

^n-4 


^n-5 

+ 

200Aco 

“n-3 

= 

^n-4 

+ 

200Aco 

^n-l 

= 

"n-3 

+ 

lOOOAco 

"n-1 

= 

^n-2 

+ 

3 

< 

o 

o 

o 

CSJ 

^n 

= 

^n-l 

+ 

3000ACO 


where j=l, 6 

where j=7, k, and k=“(n-15) 

3 

where J=k, k^, and k'>=-(k) 
where j=k 2 , and kj=2k 
where j=k 3 , .... k^, and k4=(n-4) 



APPENDIX B 


SPECIAL FEEDBACK 


The special feedback option is geared to a particular problem-the physical limit on the 
deflection angle of a rudder used in a feedback damping system. The rudder-limiting 
option, when activated, inspects the magnitude of the rudder angle. When the angle is 
less than the limit angle (RLIMIT, radians), the program functions in the 

normal fashion. At any frequency where the angle called for is greater than the limit, 
the program decreases the yaw damper gain to match the rudder’s physical limit by 
using the ratio of The phase angle remains the same as that originally 

called for. 

The option has been used to plot curves similar to curve A in the following sketch. 



Note: Stability augmentation system (SAS) 
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This information allows the engineer to determine the effect of rudder limiting on load 
exceedances. 


The special feedback option results in the following three steps in the program’s code. 

1. For each frequency, the FINDRMS subroutine SOLVEQ calculates {Q} in the same 
manner used for the standard feedback option. 

[SUMM] = [Ml] + i[gSD][Mi] + s[M 2] + ... + s^g(s)[MB] 

+ Fydamp j[Si] + s[S 2 ] + ... + s'^g(s)[S6]j 

(SUMC) = -z|f(s){C 3 } + {C 2 } 

G(Pis^ + Uis + Ri)(P2S^ + U2S + R2)(P38^ + U3S + R3KT1S + 1 ) 

Fydamp r — 

(P4S^+ U4S + R4)...(P7S^ U7S + R7KT2S + 1)(T3S + 1)...(T7S + 1 ) 

2. Now SOLVEQ solves for the rudder angle, Sr 



Sr 


NDOF 

I Fydamp (Qi(CS2iS^ + CSlis + CSOi) ) 
i = 1 


where |a| is the absolute value of the complex number a. 


3. Next, the program checks Sr against If ^R ^ ®Rum’ solution is accepted. 

If not, the feedback transfer function is factored (multiplied by 8Run,/SR) and 
steps 1 and 2 are repeated until Sr falls within the limits of 


Sim - 

However, if in five interations Sr does not converge to SRij„,> following 
message is printed and the current solution is accepted. 


*♦♦*♦♦** WARNING ERROR 16 

RUDDER LIMITING DID NOT CONVERGE ON RLIMIT( ) 
FOR FREQ( ) = 

AFTER 5 TRIES AMAG = 

RUN WILL PROCEED WITHOUT CONVERGING. 
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APPENDIX C 

MULTIPLE EXCITATION FORCING FUNCTIONS 

The analytical steps followed in section 5.0 describe the operation of FINDRMS for a 
system being excited by a single excitation source. FINDRMS can also be used to 
compute the RMS response levels and number of zero crossings for systems responding 
to multiple excitations. The multiple excitation mode, however, does not include all of 
the capabilities available during the execution of the single excitation mode. The 
following options are excluded when using multiple excitations: 

• Feedback 

• Gradual penetration 

• Built-in input spectrum 

• Kiissner lift growth function (f(s) must = 1 .) 

• Load correlation 

• Printing of loads per degrees of freedom 

• Static-elastic solution 
In addition: 

• Frequency responses and load transfer functions can be plotted for only the first 
excitation force. 

• {^2} and {C2} are null. 

The FINDRMS analytical steps for multiple excitations are: 

1. Read card input of options and constants 

2 . Read card or magnetic file input of equations of motion [Mj], [M2]>--> [^]> plus 
the excitation forces [C3]. 

Note: [C3] is a rectangular matrix with NDOF (number of degrees of freedom) 

rows and NFORC (number of excitation force) columns. The €313 element is 
the excitation scalar for the ith mode responding to the jth force. 
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3 . For each of the NFREQ frequencies, solve for [Q] the generalized coordinate 
frequency response and save [q] on "SCRATCH” file. The generalized coordinates 
are found by solving the following equations. 

[SUMM][Q] = [SUMC] 

where: 

[SUMM] = [Ml] + i[gSD][M° ] + s[M.2] + s 2 [M 3 ] 

+ g(s)[M4] + sg(s)[M5] + 8^g(s)[M6] 

where: 

s =iw = Laplace transform variable 

g(s) = Laplace transform of the Wagner function 

bjs C|S d|S 

"‘"'I's + ai s + j3| s + 7| 

[SUMC] = -z[C3] 

Note: Optionally, the input matrices [M4], [M5], [Mg], and [C3] or [^] are 

frequency dependent. Interpolation is required to define those matrices at 
each frequency, o>, before calculating [SUMM], [SUMC], and [Q]. 

Optionally, the user can delete specified degrees of freedom. This is 
accomplished by zeroing out the appropriate degrees of freedom in the 
[SUMM] and [SUMC] summation matrices, 

4 . Read ca^ or magnetic file input of load equations [M^], [M2], [Mg], and [C3]; 

where [C3] is NLD (number of load equations) by NFORC (number of excitation 
forces). 

5 . For each frequency, read [Q] (from step 3 ), calculate [SUMC], calculate the input 
spectrum [<^( 0 )^, calculate the output spectrum SPEC, optionally write data on 
files "IPLTPE” and "IRTAPE”, optionally print [SUMC], and keep a running 
integration over the frequencies of {A} and {Nq}. 

[SUMC] = [SUMM][Q] + |z|[C3] 

where: 

[SUMM] = [Ml] + s[M2] + s^[M3 ] + g(s)[M4]+ 8g(s)[M5 ] + sV8 )[M6] 

The integration for [A] and {Nq} at one load,^ , is shown as: 
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^NFREQ( SPEC,,.SPBC,, _,y^^ 


/ 1 \ fNFREQ /SPECp:n:2 + SPECp: ,2\ 

Norli^kj) N .f, ( — — — 2 >%-%!> 


where SPEC 1 1 is the output spectrum calculated for load I at frequency i and is 
found by: 


NFORC 

SPECj^i = 2 SUMCij^^(n)SUMCk^£(ni)0k,k(^i) 

k= 1 

/ NFORC NFORC ^ 

+2fi( 2 2 SUMCk^il(ni)SUMCj ji(ni)«kj(^i^ 



where {-)* denotes complex conjugate of (-), the input spectrum, and (R (-) 

denotes the real part of (•). This equation is derived from the more general form: 

NFORC NFORC 

SPECai =2 2 SUMCk.A(ni)SUMCj,i(fll)d.kj(ftl) 

k=l j=l 

by utilizing the Hermitian property of 

Note: Optionally, the input matrices [M4], [M5], [Mg], [C3] or [$] are 

frequency dependent (NKVAL > 1). Interpolation is req uired t o 
define those matrices at each frequency, a>, before calculating [SUMM], 
{SUMC},etc. 

Optionally, the user can delete specified degrees of freedom. This is 
accomplished by the zeroin g of the appropriate degrees of freedom in the 
summation matrix, [SUMM] . 

6. Print [A] and {Nq} and optionally punch them out on cards. 

7. For additional sets of load equations, repeat steps 4 through 6. 
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